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APl’LlCATlCn  Or_PATTERrJ  RECOGMITIONJPCIIMIOUPS  TO  nfiTECTIOR  OF  SEVERE  STORM 
(tYATURB  Fi^0^f  Hi:T^l)RQL0'CTC'AL  SArElLITL  DATA 

0 . Introduction 

Tho  availability  of  hicjh-resol ution  satellite  imagery  from  SMS/GOES  offers 
a unique  opportunity  to  monitor  severe  weather  phenomena  on  the  mesoscale.  The  j 

objectives  of  this  report  are  to  summarize  features  of  various  categories  of  | 

severe  storms  and  examine  the  applicability  of  pattern  recognition  techniques 
to  the  identification  and  prediction  of  severe  storms  from  visible  and  infra- 
red satellite  iniagery. 

Pattern  recognition  models  for  severe  weather  phenomena  differ  from  classi- 
cal, statistical  pattern  recognition  models  in  terms  of  feature  values,  stages 
of  decision  logic,  branching  logic,  and  feedback  capability.  Feature  values 
are  often  discrete  (representing  the  presence  or  absence  of  one  or  more  events 
associated  with  severe  weather)  rather  than  continuous.  The  decision  logic 
proceeds  through  various  stages,  often  backing  up  or  modifying  previous  deci- 
sions, until  a final  classification  (such  as  number  of  inches  of  rainfall  or 
cyclone  intensity)  is  made. 

The  degree  of  automation  which  can  be  introduced  into  these  pattern  recog- 
nition systems  can  vary  from  total  automation  to  a completely  manual  approach. 

The  cost  of  development  and/or  application  of  automatic  pattern  techniques 
must  be  considered  on  an  individual  basis  for  each  type  of  feature  to  be  ex- 
tracted. In  general,  those  features  which  must  be  obtained  by  automatically 
tracking  a sequence  of  cloud  objects  are  more  exnensive  than  those  which  can 
be  obtained  from  a single  image. 

This  investigation  began  with  a literature  survey,  a questionnaire,  on- 
site visits  to  NASA  and  NESS  and  interviews  with  meteorologists  expert  in  vari- 
ous aspects  of  severe  storm  identification  and  prediction.  Section  1 reviews 
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the  pattern  recognition  !:;oo':':  relevant  to  severe  storm  identification  and  pre- 
diction. Section  2 summarizes  for  various  categories  of  severe  weather  those 


features  wliich  citaracteri zc  each  category.  Tfie  digitized  data  necessary  to 
extract  eacfi  feature  is  considered. 

For  a completely  automatic  system  for  severe  weather  pattern  recognition 
models,  algorithms  for  (1)  cloud  classification,  (2)  cloud  segmentation,  (3) 
cloud  geometry,  and  (4)  cloud  tracking  need  to  be  developed.  A review  of  al- 
gorithms in  these  four  different  areas  and  suggestions  for  future  algorithmic 
development  are  presented  in  Section  3. 

The  concluding  section  outlines  our  recommendations  on  the  feasibility 
and  priority  of  development  of  some  pattern  recognition  techniques  for  iden- 
tification of  severe  weather  features. 
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Rt-'conni lion  Models 

One  can  consider  a number  of  different  paradicjnis  for  pattern  recognition. 

A description  of  thi?  major  models  serving  to  direct  research  on  machine  pat- 
terns appears  in  Kanal  [1974],  The  tivo  main  approaches  are  termed  (a)  the 
Feature  Extraction-Classification  model,  and  (b)  the  linguistic  or  syntactic 
model . 

Feature  Exiracti o_n  Classi ficat ion  Model 

In  the  feature  extraction  classification  model,  recognition  is  achieved 
by  making  measurements  on  the  patterns  to  be  recognized,  and  tlien  deriving 
features  from  these  measurements.  These  features  form  the  input  to  a classi- 
fication procedure  that  gives  a class,  group  or  category  assignment  for  each 
pattern.  The  available  information  from  the  pattern  environment  is  thereby 
reduced,  in  stages,  ultimately  to  a small  number  of  categories. 

A commonly  used  version  of  the  feature-extraction  classification  model 
is  one  in  which  the  'features  are  treated  as  components  of  a vector  X = 

(x]  ,X2, . . . ,Xn) . Each  pattern  is  considered  to  be  a ooint  in  the  resulting 
n-dimensional  feature  space.  Classification  is  then  treated  as  a problem  of 
dividing  this  n-dirrcnsional  space  into  exclusive  regions  Rj  , j = l,2...k, 
such  that  when  a feature  vector  falls  into  Rj , the  pattern  is  assigned  to 
class  0-  Tin’s  division  might  be  effected  on  the  basis  of  statistical  or  non- 
statistical  considerations.  Note  that  any  type  of  feature  may  be  used  in 
thi s n'odol  . 

In  the  simplest  case  where  there  is  no  variability  within  pattern  classes, 
classification  reduces  to  template  mate'nina  in  which  an  observed  pattern  is 
matched  against  a prototype  (template).  When  the  variability  is  limited,  it 
may  be  possible  to  extend  this  idea  to  the  feature  space  and  classify  patterns 

3 


according  to  tlic'ir  distance  from  the  nearest  prototype.  In  most  problems  one 
needs  to  employ  more  sophisticated  methods  of  multivariate  statistical  clas- 


sification for  defining  the  regions  R j . 


Linooistic-Syn tactic  Model 

Tfie  coiumonly  used  version  of  the  feature-extraction  classification  model 
has  been  criticized  for  focusing  primarily  on  statistical  relationships  among 
the  features  while  ignoring  other  "structural"  properties  that  seem  to  char- 
acterize patterns.  Also  the  classification  into  a region  Rj  of  n-dimensional 
space  leads  only  to  a class  designation  of  a pattern  rather  than  a description 
which  provides  some  insight  concerning  the  class,  or  which  allows  one  to  gen- 
erate patterns  belonging  to  a class. 

Various  approaches  are  being  taken  to  overcome  the  above  perceived  fail- 
ings of  the  feature-extraction  classification  model.  In  the  linguistic  model 
for  pattern  description,  patterns  are  viewed  as  sentences  in  a language  de- 
fined by  a formal  grammar.  A "primitive  extractor"  transforms  the  input  data 
into  a string  of  symbols  or  some  general  relational  structure.  The  primitive 
extractor  may  itself  be  a feature  extractor  classifier.  Then  a structural 
pattern  analyzer  uses  a formal  grammar  to  parse  the  string  and  thus  constructs 
a description  of  the  pattern. 

The  stre.ss  on  the  distinction  between  the  feature-extraction  classifi- 
cation model  and  the  linguistic-syntactic  model  misses  the  obvious  points 
that,  even  in  the  syntactic  model,  the  primitives  are  features  that  have  to 
be  extracted  from  measurements,  and  that  associating  a pattern  with  a genera- 
tive model  is  essentially  equivalent  to  classifying  the  pattern  into  categories 
represented  by  the  generative  models.  Clearly,  the  generative  models  need 
not  bo  restricted  to  formal  grammars.  One  could  include  such  familiar  models 
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as  a diffcrenlial  equation  model,  a functional  equation  model,  or  a stochas- 
tic model  such  as  a finite  state  Markov  chain  model. 

Structural  Pattern  Recogni _t;i_qn^ 

When  a formal'  model  is  not  explicitly  present,  the  terms  "ad-hoc"  or 
"heuristic"  are  used.  The  phrase  "structural  pattern  recoqnition"  refers 
to  all  pattern  recognition  approaches  based  on  defining  primitives  and  iden- 
tifying allowable  structures  in  terms  of  relationships  among  primitives,  and 
substructures  that  combine  primitives.  The  relations  might  be  boolean  expres- 
sions or  might  be  specifiable  by  some  statistical  relationships  or  by  a gen- 
erative grammar  model.  The  term  structural  pattern  recognition  represents 
less  a specific  set  of  procedures  than  an  attitude,  i.e.,  that  pattern  recog- 
nition algorithms  should  be  based  on  the  mechanisms  that  generate  and  deform 
patterns . 

Hierarchical  Classifiers 

Hierarchical  classifiers  or  decision  trees  have  been  used  extensively 
in  many  application  areas  of  multiclass  pattern  recognition.  A hierarchy 
of  classifiers  provides  a flexible  way  to  incorporate  different  feature  sub- 
sets, statistical,  linguistic,  or  ad-hoc  relationships,  and  decision  policies 
at  the  various  nodes  of  the  tree.  ; 

All  the  problems  of  feature  selection  and  classifier  design  are  accentu-  i 

ated  in  the  multiclass  multimodal  case, and  the  usual  single  stage  multivariate 
linear  or  nonlinear  discriminant  function  approach  often  fares  more  poorly 
than  anticipated. 

The  apparent  theoretical  increase  in  discriminating  ability  of  single 
stage  nonlinear  discriminant  functions  is  usually  not  realizable  in  practice  j 

because  of  problems  involved  in  estimating  the  fine  structure,  represented 
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by  tite  higher-order  relationships  among  the  features,  due  to  dimensionality 
and  sample  size  considerations  [See  Kanal  1974], 

In  an  M class  problem,  for  a given  sample  size  per  class,  dimensionality- 
sample  size  considerations  may  dictate  that  no  more  than  n features  be  used. 

It  is  likely  that  the  best  set  of  n features  to  di scrimi nante  betv/een  one 
group  of  classes  is  quite  different  from  the  best  n features  for  another  set 
of  classes.  Limiting  consideration  to  some  n features  to  be  used  in  a single 
step  to  make  an  M-v/ay  decision  may  lead  to  a forced  compromise  in  the  choice 
of  the  feature  set.  Splitting  the  decision  into  several  stages  might  be 
better.  Then  at  each  stage  of  the  decision  process,  the  n features  best 
suited  for  that  classification  task  could  be  used.  This  permits  the  multi- 
class classification  problem  to  be  decomposed  into  a series  of  less  complex 
decision  problems.  The  optimal  and  heuristic  design  of  decision  trees  v/hich 
use  both  physical  model-based  features  and  data  derived  features  is  the  sub- 
ject of  current  research  by  Dr.  Laveen  Kanal  and  his  students  at  the  Labora- 
tory for  Pattern  Analysis  of  the  University  of  Maryland  [e.g.  Kulkarni  1976, 
Kulkarni  and  Kanal  1976]. 

In  practice,  the  optimal  techniques  for  designing  decision  trees  are  not 
feasible, and  the  various  suboptimal  techniques  for  tree  design  are  best  viev.'ed 
as  subroutines  to  be  utilized  in  a recursive,  interactive  approach  to  the 
design  of  hierarchical  classifiers.  Reasons  for  taking  an  automated,  inter- 
active, graphics-oriented  approach  to  pattern  analysis  and  the  design  of 
classification  systems  were  summarized  in  Kanal  [1972],  which  also  described 
a number  of  representative  interactive  pattern  analysis  and  classification 
systems  (IPACS)  then  implemented  or  under  development. 

MIPACS,  the  Maryland  Interactive  Pattern  Analysis  and  Classification 
System  wliich  is  partially  implemented,  was  designed  to  facilitate  the  intcr- 
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active  design  of  single  stage  and  hierarchical  classifiers  [Stockman  and 
Kulkarni  (1976),  Kanal  (1977)].  MIPACS  has  been  found  very  useful  by  Ms. 

Jo  Ann  Parikh  in  interactively  designing  and  testing  decision  trees  for 
her  dissertation  on  cloud  classification  and  segmentation  [Parikh  (1977a, 
1977b)]. 

We  note  here  that  (in  unpublished  v/ork)  Scofiel'^  has  developed  a decision 
tree  for  thunderstorms  to  aid  human  decision  making  (See  Fig.]I_),  and  D. 

Tarpley  of  HESS  has  developed  a decision  tree  for  recognition  of  low  clouds. 

For  severe  storm  identification  and  recognition,  it  is  clear  chat  the 
classification  model  which  will  have  to  be  relied  u[)on  ■’•s  the  hierarchical 
classifier  or  decision  tree  with  statistical  and  other  structural  features 
and  relationships  being  used  at  the  various  nodes.  The  design  and  test- 
ing of  the  classifiers  for  some  of  the  severe  storm  features  will  have  to 
be  experimented  within  an  interactive  manner,  before  defining  an  opera- 
tional procedure  to  be  recommended  for  use  in  any  automatic  or  semi-automatic 
system. 

Some  of  the  general  problems  of  pattern  recognition,  which  affect  the 
development  of  pattern  recognition  schemes  for  a variety  of  application  areas, 
were  listed  in  Kanal  [1975].  Some  of  those  problems  impinge  upon  the  develop- 
ment of  automatic  and  semi-automatic  systems  for  severe  storm  feature  extrac- 
tion and  classification.  Additional  problems  arise  in  the  severe  storm 
area  from  the  fact  that  some  of  the  important  features  which  are  used  or 
may  be  potentially  used  to  identify  severe  storm  activity  are  not  defined  in 
a static  manner.  Rather  many  features  are  defined  by  the  nature  of  their 
evolution  in  tine,  with  the  time  involved  ranging  from  a fraction  of  an  hour 
to  days.  Procedures  for  the  automatic  extraction  of  such  time  dependent 
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features,  many  of  which  are  described  in  subsequent  sections,  are  essentially 
lacking  in  the  published  literature  on  digital  picture  processing  [Rosenfeld 
and  Kak  (1976)].  It  docs  seem  feasible  to  develop  automatic  procedures  for 
some  of  the  more  important  time  dependent  features,  and  these  are  taken  up 
in  Section  3.^,  and  in  the  recommendations.  Many  of  the  available  procedures 
for  curve  detection,  line  detection,  description  of  pattern  geometry,  and 
pattern  segmentation,  some  of  which  are  discussed  in  Section  3,  need  to  be 
examined  in  a comparative  study  for  their  relevance  to  severe  storm  feature 
detection.  The  state  of  the  art  of  image  pattern  recognition  is  unfortunately 
not  so  advanced  that  one  can  decide  without  experimentation  on  real  data 
which  procedure  is  best  for  a given  application. 
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2 . Description  of  Severe  Weather  Features^. 

A compilation  of  severe  v/eather  features  by  severe  storm  category  v/as 
prepared  from  infonration  obtained  from  a reviev/  of  the  literature,  consul- 

1 

tations  with  satellite  meteorologists,  and  responses  from  a questionnaire. 

The  format  of  the  questionnaire  is  presented  in  ttie  appendix.  The  responses 
received  are  reflected  in  the  text  of  this  report. 

Four  categories  of  severe  storms  v;ere  of  major  interest  to  the  satellite 

,! 

meteorologists:  j 

1)  Severe  windstorm  : 

2)  Heavy  rainstorm  j 

3)  Hailstorm  ^ 

4)  Tornado 

Two  other  categories  of  severe  storms--bl i zzards  and  electrical  storms  j 

(1 ightning)--wil 1 not  be  considered  in  this  study.  Blizzards  tend  to  be  | 

characterized  by  larger  scale  features  than  the  mesoscale  features  which 
have  been  used  to  identify  and/or  predict  severe  windstorms,  etc.  Further  , 

research  is  needed  before  definitive  correlations  between  specific  satellite  i 

features  and  the  occurence  of  severe  electrical  storms  can  be  made.  Many  ; 

I'l 

of  the  features  which  are  relevant  to  the  severe  storm  categories  above  j 

'j 

will  probably  prove  to  be  of  importance  in  the  case  of  electrical  storms.  j 

A list  of  severe  weather  features  I'esulting  from  merging  the  responses 

i 

to  the  questionnaire  on  this  topic  is  given  in  Table  1.  An  attempt  was  made  ? 

! 

to  place  the  most  important  features  for  each  storm  category  near  the  top  ‘ 

of  the  list. 

Most  of  tlic  features  can  bo  detected  from  observations  of  a sequence  of 
infrared  or  enhanced  infrared  images.  Included  among  these  features  arc  [ 
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splitting  cells,  merging  cells,  speed-of-cel 1 movement,  vertical  growth 
rate  of  colls,  change  in  top  heights  and  life  cycles  of  overshooting  tops. 

Wind  fields  which  can  be  automatically  extracted  from  visible  or  infrared  se- 
euencos  by  cloud  tracking  algorithms  can  be  used  to  determine  low-level  vor- 
ticity,  patterns  of  upper-level  divergence  combined  with  low-level  conver- 
gence, whether  a storm  is  deviating  to  the  right  or  left  of  the  mean  v/ind, 
and  upper-level  jets  in  the  vicinity  of  storms.  Specific  cloud  types  or 
configurations  associated  with  these  features  also  can  be  used  to  identify 
their  presence  in  a single  imago.  Mesohigh  boundary  interactions,  pendant- 
shaped cells,  low-level  cumulus  cloud  line  feeder  bands,  the  texture  of  a 
cell,  and  the  location  of  the  steepest  temperature  gradient  are  severe  weather 
features  which  can  be  observed  in  a single  image.  Overshooting  tops  are,  in 
general,  too  small  to  be  observed  in  the  infrared  images  where  the  resolu- 
tion is  4 nautical  miles  (n.m.)  hut  can  be  detected  in  the  visible  images 
whore  the  t'esolution  is  1/2  n.m.  or  1 n.m. 

A brief  description  of  eacn  of  the  features  in  Table  1 is  included  in 
tliis  section.  Algorithms  to  extract  each  feature  arc  considered 
in  Section  3.  Section  3 provides  insight  into  the  na- 

ture of  the  problems  which  must  be  solved  if  a completely  automatic  system 
IS  U)  be  developed. 

2.1.  Re s 0 h i g h_ F. o_^i i d a ry  Interact i_o n 

The  intersection  or  merger  of  the  outer  boundary  of  a mesoscale  hign 
pressure  system  with  another  convective  boundary  (such  as  a front,  squall 
line,  or  another  mesohigh)  was  shown  by  Purdom  [74,76a]  to  almost  always  re- 
sult in  an  increase  in  convective  activity.  Rain,  a decrease  in  tempera- 
ture, a pressure  surge,  gust  fronts,  and  tornado  activity  may  occur  at 


10 


I 

I 

these  points  ot  convective  intersection.  There  are  two  basic  cloud  pat- 
terns described  by  Gurka  [76]  which  are  associated  with  mesohighs.  In  the 
first  case,  the  outer  boundary  of  the  mesohigh  appears  "as  an  oval  shaped 
leading  edge  to  the  thunderstorm  wi th  the  anvil  cirrus  curving  anticycloni- 
cally  out  of  the  rear  of  the  (thunderstorm)  cluster".  In  the  second  case, 
the  outer  boundary  appears  as  "an  arc  shaped  line  of  convective  clouds  ad- 
vancing outward  from  a dissipating  thunderstorm  area".  The  parent  cumulonim- 
bus cell  may  exhibit  a "comma"  tail. 

The  comma-shaped  cloud  pattern,  which  is  associated  with  areas  of  upper 
tropospheri G.  vorti ci ty , is  often  indicative  of  severe  weather.  Mathews  and 
Johnston  [76]  reported  the  occurrence  of  tornado  activity,  gust  storms  in 
excess  of  50  knots,  and  hail  in  the  area  of  a compact,  wel 1 -organized  comma 
cloud  system.  They  attributed  the  severity  of  the  thunderstorm  outbreak  to  the 
relationshio  between  well-defined  comma  cloud  systems  and  concentrations  of 
positive  vorticity  advection  (strong  upward  vertical  motion).  Parmenter  [76] 
described  a rain/snow  producing  sub-synoptic  system  in  which  the  cloudiness 
associated  with  a surface  wave  had  a distinct  comma  shape  with  colder  tops 
along  the  northern  side. 

In  order  to  determine  or  differentiate  areas  of  severe  weather  along 
the  mesohigh  boundary,  Gurka  [76]  recommended  in  the  case  of  gust  fronts  use 
of  enhanced  infrared  and  visible  imagery  to  locate  features  such  as  tempera- 
ture cf  the  coldest  tops  in  tiie  thunderstorm  cluster  behind  the  mesohigh 
boundary,  the  presence  of  overshooting  tops,  cloud  edge  gradients,  and  speed 
of  motion  of  the  arc  bounding  the  mesohigh.  The  role  of  intersecting  arc 
linos  (mesohigh  boundary  interactions)  was  also  emphasized. 
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2.2.  Level  Cumulus  C l_o  ir!_  M n e_  JT  e e d e r_  I Li  n d s 

Low-level  cumulus  cloud  line  feeder  bands  refer  to  the  mesoscale  lines 
of  small  cumulus  clouds  which  laove  or  feed  into  a thunderstorm  system,  main- 
taining and  replenishing  the  moisture  of  the  thunderstorm  system.  These 
cloud  lines  represent  areas  of  low-level  moisture  convergence.  When  these 
convective  lines  intersect  or  merge  into  another  convective  boundary  or  lino, 
enhanced  thunderstorm  activity  almost  always  results.  An  example  of  this 
type  of  merger  appeared  in  Purdom  [74].  The  location  of  now  thun.dcrstorms 
along  the  cloud  line  of  an  old  nosohigh  boundary  was  sliown  to  coincide  with 
the  point  where  low-level  cumulus  lines  merged  into  it  from  the  south. 

The  recognition  of  convective  lines  in  satellite  imagery  plays  a major 
role  in  the  prediction  of  thunderstorm  activity.  Oliver  and  Purdom  [74] 
stated  that  "almost  all  new  large  convective  activity  in  tiie  tropics  forms  on 
a pre-existing  line".  They  differentiated  various  types  of  line  patterns 
in  satellite  imagery  including  synoptic  frontal  bonds,  squall  lines,  rope- 
like clouds,  and  mesoscale  lines  of  small  cumulus  clouds.  Convective  lines 
representing  areas  of  low-level  moisture  convergence  may  be  terrain  related. 
However,  Purdom  [74]  noted:  "Regardless  of  the  generation  mechanism  for  a 

convective  boundary  or  line,  when  it  merges  with  and  intersects  another  con- 
vective boundary  or  line,  more  intense  convective  activity  alinost  always  re- 
sults. " 

2.3.  Steepest  T emnerature  Gra cHjni t 

The  location  of  a gust  front  along  the  leading  edge  of  a cumulonimbus 
system  and  the  area  of  most  significant  precipitation  in  a convective  syste;'- 
can  be  determined  by  analysis  of  the  temperature  gradients  witiiin  the  thun- 
derstorm systems.  Gurka  [76]  noted  that  "tiie  gust  front  is  usually  located 
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very  close  to  the  strongest  temperature  gradient  in  the  infrared  pictures, 
near  the  leading  edge  of  the  cumulonimbus".  Oliver  and  Scofield  [76]  ob- 
served that  the  area  of  most  significant  precipitation  (which  may  consist 
of  only  one- tenth  or  less  of  the  entire  anvil  area)  in  a convective  system 
occurred  in  the  upwind  [)ortion  of  the  system.  The  upwind  portion  of  the  con- 
vective system  can  be  determined  in  a single  enhanced  infrared  image  from  the 
shape  of  the  individual  anvils  and  t!ic  location  of  the  steepest  temperature 
gradient.  Anvil  cirrus  has  a distinct  sharp  edge  on  the  upwind  side  and  a 
fuzzy  edge  downwind.  The  plumes  often  extend  downwind  several  degrees. 

The  location  of  the  largest  temperature  gradients  in  images  formed  by  sum- 
ming digital  values  for  enhanced  infrared  images  over  a 6-hour  period  can 
also  be  used  to  determine  the  upwind  portion  of  the  convective  system.  A 
high  degree  of  correlation  (correlation  factor  of  0.88)  was  obtained  by  Oliver 
and  Scofield  [76]  between  observed  6-hour  precipitation  and  cumulative  digi- 
tal infrared  values  in  the  upwind  portion  of  the  convective  system. 

Enhanced  infrared  imagery  is  used  to  "better  display  low  clouds,  to 
better  identify  the  structure  of  convective  storms  and  to  locate  sea-surface 
temperature  gradients"  (Anderson[74] ) . Enhanced  infrared  images  are  ob- 
tained by  transforming  the  gray-scale  of  infrared  images  by  application  of 
specific  eniiancemcnt  curves.  Display  devices  for  infrared  imagery  are,  in 
general,  limited  from  16  gray  shades  (or  less  for  facsimile  recorders)  to  64 
gray  shades.  If  infrared  imagery  is  not  enhanced,  temperature  increments  are 
equal  for  each  shade  of  gray  in  the  display.  In  order  to  enhance  tlie  data, 
the  range  of  possible  temperature  values  is  partitioned  into  temperature  in- 
tervals or  segments  and  for  each  segment  a linear  mapping  from  the  tempera- 
ture values  witliin  the  segment  to  a specific  gray  scale  range  is  specified. 
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If  no  enhancement  is  desired  for  the  segment,  the  linear  function  consists 
of  a straight  line  with  slope  factor  1.  The  enhancement  curve  may  assign 
gray  shades  out  of  sequence  instead  of  varying  gray  shades  from  black  to 
white  as  the  temperature  decreases.  This  type  of  enliancement  (alternating 
light  and  dark  shades)  is  particularly  effective  for  identifying  the  struc- 
ture of  convective  storiDs.  The  enhanced  infrared  pictures  used  by  Gurka  [7'1] 
were  produced  by  assigning  gray  shade  values  to  temperature  intervals  as 
specified  below: 


+30°C  to  -44°C  --  linear  assignment  from  black  to  near  white 

-45^C  to  -57°C  — dark  gray 

-58°C  to  -63'’C  --  light  gray 

colder  than- -63°C  --  black 

The  following  enhancement  procedure  was  used  by  Oliver  and  Scofield  [76]  for 
sum.mer  convection  in  the  middle  latitudes: 


+17°  to  -3°C 
-3°  to  -23°C 
-24°C  to  -43°C 
-44°C  to  -5S°C 
-59°C  to  -63°C 
-64°C  to  -68°C 
below  -68°C 


— dark  gray 
--  light  gray 

— white 

— dark  gray 

— light  gray 
--  blackest 

— whitest 


A comparison  of  three  eniiancemoiit  curves  for  GOES-1  infrared  imagery  avail- 
able in  the  k'ashi ngton  Satellite  Field  Services  Station  on  February  1976  can 
be  found  in  Clark  [’6]. 

The  location  of  the  steepest  temperature  gradient  in  enhanced  infrared 
imagery  is  determined  by  observing  a narrowing  of  the  width  or  distance  be- 
tween gray  shade  contours.  According  to  Scofield  [76],  "There  arc  situations 
when  the  upwind  portion  (area  of  tight  IR  temperature  gradient)  and  active 
thunderstorm  clusters  are  difficult  to  find  because:  (1)  the  tight  IR  temp- 

erature gradients  cannot  be  discerned,  or  (2)  the  colder  enhanced  IR  contours 
overspread  a large  area."  There  is  clearly  a need  to  incorporate  pattern  re- 
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cognition  tocftniques  for  analysis  of  inirarcd  gradients  into  a sys^o'a  for 
producing  enhanced  infrared  imagery.  The  need  to  develop  pattern  recogni- 
tion techniejues  that  enhance  and  identify  the  surface  thermal  gradients  v/as 
discussed  by  Miller  et  al.  [7h]  under  the  heading  of  "Severe  Storm  Potential 
Analysis".  They  recommended  use  of  thermal  gradient  maps  to  quantify  phen- 
omena such  as  the  motion  of  arc  clouds  and  to  delineate  "tlie  relative  hot 
and  cold  'tongues'  that  seem  to  be  necessary  conditions  for  certain  kinds 
of  severe  weather  generation". 

2.4.  '[opHej^hts^ 

"The  highest  and  coldest  clouds  form  where  the  thunderstorms  are  most  vig- 
orous" (Oliver  A Scofield  [76]).  Near  the  convective  cells  with  the  coldest  tops 
tlie  strongest  winds,  the  heaviest  rain,  and  the  most  severe  thunderstorm 
activity  occur.  An  outbreak  of  severe  weather  in  the  Midwest  which  pro- 
duced 4 tornadoes,  14  occurrences  of  gusts  in  excess  of  50  knots,  wind  damage, 
and  hail  up  to  coif  ball  size  was  found  by  Mathews  and  Johnston  (76]  to  be 
associated  with  a comma  cloud  pattern  in  which  cloud  tops  up  to  46,000  ft. 
were  indicated  by  radar. 

The  height  of  cloud  tops  can  be  obtained  f roai  the  cloud-top  temperatures 
and  a vci'tical  temperature  profile  of  the  area  of  interest.  Procedures  to 
determine  cloud-top  temperatures  from  satellite  radiance  measurements  gen- 
erally make  two  basic  assumptions:  (1)  clouds  arc  opaque  and  hence  their 
emissivi tics  are  equal  to  unity,  and  (2)  a constant  correction  factor  may  be 
applied  for  water  vapor  contributions  above  the  cloud.  Liou  [75  ] noted 
that  both  assumptions  seem  justified  for  cumulonimbus  clouds  which  are  fairly 
high  in  the  atmosphere  but  added  that  "it  seems  unlikely  that  a cloud  whose 
emissivity  equals  unity  may  be  isolated  from  satellite  observations",  and 
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that  "unless  a reliable  radiative  transfer  calculation  can  be  carried  out, 
the  uncertainties  of  the  einissivity  of  thick  clouds  or  multilayered  clouds 
may  lead  to  serious  errors  in  the  evaluation  of  cloud-top  temperatures."  If 
the  field  of  view  of  the  satellite  sensor  is  not  completely  filled  by  clouds 
which  arc  sufficiently  dense  to  shield  the  satellite  sensor  from  radiation 
below  the  cloud,  the  radiance  at  the  satellite  will  represent  the  sum  of 
radiance  from  the  cloud  and  radiance  from  surfaces  beneath  the  cloud.  If 
other  types  of  clouds  are  also  in  the  field  of  view,  for  example,  layer-type 
clouds  such  as  cirrostratus  and  altostratus,  the  assumption  of  an  emissivity 
of  1.0  may  not  be  valid.  Discussions  of  relationships  between  cloud-top 
temperatures,  brightness,  and  heights  can  be  found  in  Mosher  [76],  Shenk 
and  Curran  [73],  Park  et  al.  [74],  Gruber  [75],  Liou  [75],  and  Young  [75]. 

The  cloud  height  program  of  Mosher  [76]  determines  the  height  of  wind  tracer 
clouds  by  using  visible  brightness,  solar  zenith  angle,  satellite  zenith 
angle,  relative  azimuth  angle,  cloud  type  (ice  or  water  clouds),  etc.,  to 
calculate  cloud  emissivity  as  a function  of  the  optical  thickness  of  a cloud, 
further  investigation  is  needed  to  determine  whether  or  not  within  a thunder- 
storm system  areas  of  coldest  satellite  radiance  measurements  can  be  equated 
with  areas  of  coldest  cloud  tops  (i.e.,if  the  assumption  of  cloud  emissivity 
equal  to  unity  suffices)  or  if  more  complicated  algorithms,  such  as  the  cloud 
height  algoritlim  of  Mosher  [76],  need  to  be  implemented. 

2.5.  Vertical  Grow th  Rate  of  Cell 

The  vertical  growth  rate  of  a cell  is  determined  from  enhanced  infrared 
imagery  by  observing  the  change  of  satellite  radiance  measurement  with  time. 
Assuming  that  the  changes  in  radiance  measurements  parallel  changes  in  tem- 
perature and  that  temperature  decreases  with  height,  then  the  expansion  of  a 
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gray  shade  contour  representing  the  coldest  radiance  measurements  denotes 
a vertical  growth  of  the  cell.  The  rate  of  growth  is  represented  l’)y  the 
change  in  area  of  this  contour  over  ttio  time  period  of  interest.  This  change 
will  be  a function  of  t.he  particular  enhancement  curve  used  to  produce  thC' 
enhanced  iiiagcry. 

The  amount  of  expansion  (measured  in  terns  of  degrees  of  latitude)  of 
cloud  tops  in  t!ie  up'vind  portion  of  a convective  system  was  one  of  the  fac- 
tors in  Scofield's  decision  tree  nrocedure  (Fig.  I)  which  contributed  toward  an 
increase  in  the  amount  of  convective  rainfall.  Larger  amounts  of  pi'ecipi- 
tation  ware  associated  witii  colder  tops  and  with  rapidly  expanding  tops. 

Winds  of  50  knots  were  reported  in  the  area  of  a cloud  pattern  typical  of 
air  mass  thunderstorms  when  the  parent  cumulonimbus  cell  was  rapidly  grow- 
ing and  producing  heevy  rainfall  (Gurka  L^-T]). 

2.6.  _0y^s_h 0 o tine  Too s 

Overshooting  tops  refer  to  thunderstorm  tops  which  "shoot  over"  or 

k 

above  the  arvil  portion  of  the  IhLinderstorm.  The  upper  limtt  of  thunderstorm 
growth  is  determined  by  the  height  of  the  stratosphere . The  tropopause, 
the  base  of  the  stratosphere,  is  defined  as  "the  level  at  which  temperature 
begins  to  increase  vritii  hcigiit  cr  at  least  decreases  at  a rate  below  1.1"F. 
per  1000  ft.  (r.at.L.'n  [61]).  '.."non  the  cloud  air  moves  into  the  stratosphere, 

it  \.’ill  be  subjected  to  a downward  force  since  the  ascending  air  is  colder 
and  hence  more  dense  than  the  surrounding  air.  Tlie  updraft  speeds  in  the 
thunderstorm,  however,  may  be  so  high  that  the  upward  r;omentum  will  propel 
the  cloud  air  several  thousand  feet  above  tiic  tropooause  until  the  downward 
force  brings  it  to  a lualt. 


Overshooting  tops  can  be  recognized  fron  both  visible  and  infrared 
imagery.  The  presence  of  overshooting  tops  in  infrared  imagery  can  be  as- 
certained by  determining  if  the  cloud-top  temperature  is  higher  than  the 
temperature  of  the  tronospause  at  the  given  location.  Killer  et  al  [75] 
suggested  that  overshooting  tops  may  be  "delineated  by  simple  gradient  or 
curvature  enhancement  techniques".  A second  indication  of  overshooting 
tops  is  the  appearance  of  an  uneven  or  lumpy  texture  vjithin  the  thunderstorm 
system.  Individual  overshooting  tops  are  usually  too  small  to  be  seen  in 
infrared  imagery  (4-mile  resolution).  If  visible  imagery  (1/2-  or  1- 
milc  resolution)  is  available,  overshooting  tops  can  be  detected  by  the  shad- 
ow they  cast  upon  the  anvil. 

Overshooting  tops  are  associated  with  the  most  hazardous  parts  of  thun- 
derstorms. If  a cloud  is  capable  of  producing  large  hailstones,  it  must 
have  very  strong  updrafts.  These  updrafts  would  furnish  the  upward  momen- 
tum to  produce  overshooting  tops.  According  to  Battan  [61  ] , if  a hailstone 
is  to  reach  a diameter  of  3 inches,  it  must  take  several  trips  to  the  top 
of  a tall  thunderstorm  and  down  again.  Since  a 3-inch  hailstone  has  a fall 
speed  of  about  6000  feet  per  minute,  vortical  motions  of  this  magnitude  must 
exist  within  the  cloud.  When  analyzing  cloud  patterns  associated  with  strong 
wind  zones,  Gurka  [74  ] stated  that  "on  visible  imagery,  the  strongest  cells 
can  sometimes  be  located  by  overshooting  tops  above  the  anvil."  Scofield's 
decision  tree  procedure (Fig. 1 ) for  half  hourly  estimation  of  convective  rain- 
fall increases  the  half  hourly  estimated  amount  of  precipitation  by  0.50 
inch  if  overshooting  tops  (located  on  visible  imagery)  occurred  in  the  up- 
wind portion  of  a thunderstorm  system.  For  an  example  of  visible  imagery 
showing  a line  of  overshooting  tops  in  a thunderstorm  cluster,  see  Scofield 
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2.7.  i c h I Textu r cd 

Thundorstorm  clmul'j  v.'hich  produco  the  hcdviost  precipitation  appear 
bripht  and  textured  in  ttie  visible  imagery,  bright,  smooth  clouds  are  of- 
ten representative  of  middle  and  high  cloud  debris  v/ithoul  precipitation 
(Src'-field  [76]).  Darker  clouds  and  translucent  clouds  also  are  not  generally 
indicative  of  heavy  amounts  of  falling  precipitation.  In  addition  to  a 
bri ght-textured  appearance,  features  which  can  be  extracted  from  visible 
imagery  for  characterization  of  heavy  rainstorms  include  overshooting  tops, 
merging  tfiundcrstorms  (also  can  be  observed  in  enhanced  infrared),  and  low- 
level  cumulus  cloud-line  feeder  bands. 

2.8.  Upp  e r^  to  v cj_  1 c j n vt  y__o  f _^qrn 

A jet  stream  is  defined  by  the  liorld  Meteorological  Organization  (K'MO) 
as  "a  strong  narrow  current,  concentrated  along  a quasi-horizontal  axis  in 
the  upper  troposphere  or  in  the  stratosphere,  characterized  by  strong  verti- 
cal and  lateral  wind  shears  and  featuring  one  or  more  velocity  maxima". 

There  are  two  main  westerly  jet  streams  at  upper  troposphere  level:  (1) 

the  subtropical  jot  stream  which  is  found  betv;een  30,000  and  40,000  feet 
above  sea  level  between  30”  and  35°  latitude  and  (2)  the  circumpolar  or  polar 
jet  stream  located  somewhat  lower  in  the  atmosphere  and  further  away  from 
the  equator  than  tlie  subtropical  jet  stream.  According  to  the  definition 
adopted  by  the  Viorld  Meteorological  Organization;  "Normally  a jet  stream  is 
thousands  of  kilometers  in  length,  hundreds  of  kilometers  in  width  and  some 
kilometers  in  depth.  The  vertical  shear  of  wind  is  of  the  order  of  5 to  10 
mps  (1  meter  per  second  - 2 knots)  per  kilometer,  and  the  lateral  shear  is  of 
the  order  of  5 mos  per  100  km.  An  arbitrary  lower  limit  of  30  mps  is  as- 
signed to  the  speed  of  the  wind  along  the  axis  of  a jet  stream".  A jet 


axis  is  the  axis  of  strongest  v;inds  within  a jet  stream.  Jet  maxima  (velo- 
city maxima)  are  centers  of  liigii  wind  speed  .ilong  the  jet  streams.  Winds 
may  reach  more  than  250  knots  along  a jet  stt'eam  axis. 

Location  of  a jet  axis  can  I'C  determined  from  satellite  derived  upper- 
level  wind  data  or  from  examination  of  cloud  i>al!erns  in  the  vicinity  of 
the  jet  stream.  Auvinc  and  Sikda;  [73]  oL'.erved  that  "the  jet  axis  closely 
follows  the  zero  vorticity  line*  with  ant i; '/•  l-amc  vorticity  to  tlie  southeast 
and  cyclonic  to  the  northwest",  whiff  "was  ir;  ,i(-.'rd  with  the  theoretical 
expectation  for  a westerly  jet  streai,;".  They  u-.ed  M mruso  and  Lndlich's  [73] 
Analysis  Program  to  obtain  a comouter  analysis  of  the  diverg- 
ence and  relative  vorticity  fields.  "The  main  jet-stream  cloud  features  are 
long  shadow  lines,  large  cirrus  shields  with  sharp  boundaries,  long  cirrus 
bands,  cirrus  streaks,  and  transverse  bands  witliin  cirrus  cloud  formations" 
(Anderson  [7^]).  In  order  to  locate  on  visible  imagery  the  [)olar  jet  (which 
breeds  blizzards  and  stormy  weather  over  the  United  States),  one  should  look 
for  the  sharp  poleward  edge  of  a largo,  slightly  anticyclonically-curved 
cirrus  shield.  The  shadow  cast  by  the  high  jet-stream  cirrus  clouds  on 
lower  clouds  or  surface  features  usually  extends  far  from  the  edge  of  the 
cloud.  "Poleward  from  the  jet  streams  the  clouds  present  a 'lumpy'  appear- 
ance, while  on  the  equatorward  side  the  clouds  look  much  more  even"  (Ander- 
son [74]).  Relationships  betv/een  synoptic  features  (such  as  jet  streams) 
and  representative  cloud  patterns,  although  generally  valid,  do  not  always 
apply  in  each  specific  instance.  For  example,  Doswell  and  Schaeffer  [76]  re- 
ported a finding  on  May  4,  1976,  in  which  the  position  of  cirriform  bands  and 
streaks  was  well  south  of  the  jet  core  which  was  "at  variance  with  the  litera- 
ture on  satellite  interpretation." 


Middle  and  upper-level  jets  are  often  associated  with  outbreaks  of 
severe  weather.  Anderson  [ 7A]  pointed  out  that  "an  important  factor  in 
forecastinq  severe  weather  is  the  presence  of  a mechanism  that  produces  upper- 
level  divergence,"  and  "in  most  cases  this  is  a jet  stream."  Tornado 
outbreaks  have  occurred  near  the  intersection  of  a low-level  and  upper- 
level  jot  (Miller  [71])  and  near  the  intersection  of  a dry  line  and  a jet 
stream  (Miller  [76]).  According  to  Anderson  [7A],  "tornadoes  occurring  in 
the  presence  of  the  subtropical  jet  stream  are  normally  beneath  its  inter- 
section with  the  squall  line  and  to  the  north  of  the  jot  axis."  Scofield 
[75]  traced  a fioavy  snowstorm  to  an  unper  air  disturbance  in  the  vicinity 
of  a 300-n'b  jet  axis.  Auvine  and  Sikdar  [73]  located  the  jet  core  to  the 
west  of  storm  regions  which  produced  severe  weather  including  tornadoes  and 
located  the  jet  axis  to  the  northwest  of  the  nreci pi taticn  area.  In  the 
cases  they  studied,  they  found  similar  fluctuations  between  cirrus  tracer 
r’otions  near  the  jet  stream  and  nearby  severe  storm  intensity  which  suggested 
"some  sort  of  direct  interactions  of  the  severe  storm  evolutions  with  the 
jet- stream  flow  field." 

2.9.  Pendant- Shaped  Cel  1 s 

Pendant-shaped  cells  or  "carrot"  cells  refer  to  thunderstorm  cells  which 
develop  in  a strong, verti cal  shear  environment.  These  thunderstorm  cells, 
which  arc  shaped  by  the  veer  of  wind  with  height,  are  stretched  out  horizon- 
tally so  that  the  elongated  axis  of  the  anvil  is  parallel  to  the  vertical 
wind  shear  between  the  lower  and  upper  troposphere.  The  vertical  wind  shear 
vector  is  parallel  to  the  mid-troposplieri c isotherms.  Pendant-shaped  cells 
often  occur  in  close  proximity  to  jet  streams;  typically,  in  the  vicinity  of 
a loiv-level  southerly  jet  and  a high-level  westerly  jet.  The  wind  shift 
from  southeast  to  south  to  southwest  with  height,  produces  the  characteri Stic 
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pendant  shape  of  the  cells. 

Anderson  [7^]  pointed  out  that  "veering  of  v/inds  v/ith  height  is  an  im- 
portant parameter  in  forecasting  severe  v;eathcr."  NOAA-1  infrared  imagery, 
depicting  pendant-shaped  cells  in  an  area  in  which  tornadoes  occurred  be- 
neath the  intersection  of  low- and'mid-level  jots,  can  be  found  in  Figure  5- 
E-3  of  Anderson  [74],  Although  pendant-shaped  colls  are  generally  associated 
with  thunderstorms,  their  significance  in  determining  the  severity  of  a 
thunderstorm  is  questionable.  Further  statistical  analysis  is  needed  to 
evaluate  the  importance  of  including  pendant-shaped  cells  as  a feature  in 
severe  storm  pattern  recognition  models. 

2.10.  1^' qht  (o r Left)  Deviating  Storms 

Right  (or  loft)  deviating  storms  are  thunderstorm  cells  which  move  to 
the  right  (or  left)  of  the  mean  tropospheric  wind.  Right  moving  (right 
deviating)  storms  often  occur  after  a thunderstorm  cell  splits  into  two  sep- 
arate sections,  with  one  section  (left  deviating)  moving  to  the  left  of  the 
mean  wind  and  one  section  (right  deviating)  to  the  right.  The  right 
deviating,  cyclonically  rotating  cell  is  associated  with  severe  thunderstorms. 
"One  of  the  characteristics  of  a convective  system  is  for  thunderstorm  gen- 
eration to  take  place  on  the  upwind  side  and  to  propagate  to  the  right  of 
the  mean  tropospheric  wind"  (Scofield  [70]). 

The  movement  or  deviation  of  a storm  from  the  mean  wind  can  be  observed 
from  sequences  of  enhanced  infrared  imagery.  The  coldest  temperatures  of 
the  right  deviating  thunderstorm  move  to  the  right  of  the  mean  motion.  The 
moan  wind  or  wind  direction  can  be  determined  by  tracking  appropriate  cloud 
targets  for  wind  velocity  estimation.  If  the  cliange  in  wind  direction  with 
height  is  minimal,  the  wind  direction  will  be  parallel  to  and  can  be  esti- 
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mated  from  the  orientation  of  the  cirrus  blowoff  (cirrus  plumes).  However, 
this  proccdui'e  should  not  be  used  in  an  environment  of  strong  vertical  wind 
shear. 

In  an  anal y'’- is  ot  an  outbreak  of  severe  weather  over  the  midwestern 
United  States  on  May  5,  1971,  Anderson  [74]  pointed  out  that  "the  threat 

tracks  of  the  tornado  ''roducing  storms  are  all  oriented  at  a far  greater 

angle  to  the  right  of  the  low-level  flow  than  those  of  the  non-tornado  pro- 
ducers", which  "indicates  that  the  tornado-producing  thunderstorms  were  all 
cyclonically  rotating,  right  deviating  severe  storms."  Burgess  and  Brown 
[73]'  c'f'i  analysis  of  Doppler  radar  data  from  a right  moving  thunderstorm, 
observed  that  the  right  '...ving  cell  (referred  to  as  the  Davis  storm)  pos- 
sessed all  of  tlie  radar  supercell  characteristics.  "At  least  five  short- 

lived tornadoes  occurred"  in  the  Davis  storm  area. 


2.11.  ]^w-_Lev_el^  Vortjpcjity^ 

Positive  low-level  vorticity  and  mid-level  positive  vorticity  advection 
are  indicative  of  severe  weather  which  produces  tornadoes.  Vorticity  is  a 
vector  quantity  V.  which  measures  rotation.  If  Hy,  and  sre  the 

components  of  f2  and  u,v,  and  w are  the  components  of  the  velocity  vector  V 
with  respect  to  axes  x.y,z,  which  are  fixed  in  the  earth.,  then  the  relative 
vorticity  Qz  is  defitied  as  below; 
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The  term  vorticity  (or  relative  vorticity)  usually  refers  only  to  the  vorti- 
cal component  which  measures  spin  in  the  horizontal  plane.  A parcel  of 
air  which  rotates  in  the  same  sense  as  the  Earth  does  in  space  is  said  to 
have  positive,  or  cyclonic,  vorticity.  In  order  to  comoute  vorticity  "all 
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that  need  be  measured  is  the  velocity  of  the  air.  If  this  could  be  done  with 
sufficient  accuracy  over  a larqe  area,  it  v/ould  be  possible  immediately  to 
construct  charts  showing  how  the  vorticity  of  the  wind  is  distributed" 

(Sutton  [61]). 

low-level  cumulus  clouds  which  exhibit  cyclonic  sfioar  and  curvature  are 
indicative  of  positive  low-level  vorticity.  Positive  vorticity  advection 
"can  be  identified  where  there  is  a small  area  of  enhanced  convective  activ- 
ity or,  in  certain  cases,  a formation  in  the  shape  of  a comma"  (Anderson  [74]). 
"Cellular  cloud  patterns  aid  in  the  identification  of.. regions  of  positive 
vorticity  advection."  According  to  Clark  [76]  positive  vorticity  advccting 
over  an  area  at  mid-tropospheric  levels  can  be  "best  defined  on  IR  imagery 
as  middle-level  alto  clouds  as  opposed  to  higher,  colder  cirrus  clouds." 

Positive  low-level  vorticity  has  been  shown  by  Charba  [76]  to  be  corre- 
lated with  the  outbreak  of  tornadoes.  Parnenter  [76]  observed  from  movie 
loops  of  visible  SMS-2  imagery  a low-level  circulation  with  the  comma  pat- 
tern associated  with  rain/snow  producing  sub-synoptic  systems  around  a large 
Gulf  of  Alaska  low.  The  compact  and  well-defined  nature  of  a vorticity  com- 
ma cloud  system,  observed  by  Mathews  and  Johnston  [76]  was  believed  by  the 
authors  to  imply  "a  concentration  of  positive  vorticity  advection  (strong 
upward  vertical  motion),  and  thus  [be] . . . indi cati ve  of  the  severity  of  this 
thunderstorm  outbreak."  The  importance  of  mid-tropospheric  vorticity  ad- 
vection in  indicating  waves  on  fronts  which  were  associated  with  areas  of 


I heavy  precipitation  was  illustrated  in  an  analysis  of  frontal  cloud  systems 

by  Parmenter  and  Anderson  [74]. 
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2.12.  lj'[?|jor-Lpvo  1 Divcroor.cp  and  l.ovj-Lovcl  Convornonce 

ConvcM'Qcnco  and  divergence  refer  to  the  accumulation  and  depletion 
(respectively)  of  mass  in  a volume  of  fluid.  Low-level  horizontal  conver- 
gence implies  that  air  is  crowded  into  a smaller  horizontal  area  and,  since 
it  cannot  accuniulate  there,  is  forced  to  move  upwards.  As  the  air  ascends, 
it  cools  adi abatical ly  and  condenses  its  water  vapor  into  clouds,  rain,  and 
snow.  E.xcept  in  small-scale  disturbances  such  as  tornadoes,  the  rate  of 
change  of  pressure  at  the  surface  is  small  ("barometric  tendency")  whicli 
means  that  the  total  divergence  in  a column  of  air  must  be  small.  Therefore, 
in  order  to  "balance  a large  convergence  of  air  in  the  lower  levels  of  a 
cyclone,  there  must  be  a slightly  greater  divergence  aloft"  (.Sutton  [61]). 

The  matliematical  term  "divergence"  includes  both  the  meteorological  con- 
cepts of  convergence  and  divergence.  If  u,v,  and  w are  the  components  of 
the  velocity  vector  V in  a coordinate  system  with  a.xes  x,y,  and  z,  then  the 
horizontal  divergence  of  V is  defined  as 
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For  positive  diV[|  V,  divergence  is  said  to  occur  and  for  negative  diV||  V,  con- 
vergence. 

The  presence  of  convergence  or  divergence  is,  in  practice,  often  esti- 
mated from  the  confluence  or  diffluence  of  the  streamlines.  Areas  of  con- 
vergence may  appear  in  satellite  imagery  as  very  fine,  closely  spaced  cloud 
bands  (r'olroyd  [71]).  Interoretations  of  cellular  cloud  patterns  may  enable 
the  meteorologist  to  identify  divergent  flow  in  the  cold  air  behind  polar 
fronts  over  o''oanic  areas  (Anderson  [7^]). 

Auvine  and  Sikdar  [7.7]  pointed  out  that  "divergence  in  the  upper  atmo- 
sphere gives  a general  indication  of  the  region  of  upward  motion."  They  con- 
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eluded  tiuit  "fields  of  divergence  and  relative  vorticity  obtained  from  the 
(wind)  vectors  exhibit  a consistent  pattern  with  reference  to  the  location 
of  the  jet  axis  and  thunderstorm  prec i i)i taticn  area".  They  noticed,  however, 
that  although  the  precipitation  area  was  associated  v/ith  an  area  of  diver- 
gence, the  location  of  maximum  divergence  did  not  coincide  exactly  with  the 
location  of  the  severe  storiii  precipitation  area.  Among  the  features  ob- 
served by  Scofield  p5]  when  tracking  an  upper  ai r disturbance  which  produced 
heavy  snowstorms  were  I'-ositive  vorticity  advection  and  upper  air  divergence. 

Storms  are  often  reported  in  the  area  of  convergence  of  land,  ocean,  and 
sea  breezes.  A discussion  o""  a variety  of  factors  which  influence  the  de- 
velopment of  the  sea  breeze  can  be  found  in  Purdom  [74].  Local  areas  of  con- 
vergence of  tiio  sea  breeze  lead  to  a strengthening  of  the  cumulus  activity 
along  the  sea  breeze  front  with  thunderstorms  often  forming  along  the  merging 
of  convective  lines  such  as  those  resulting  from  two  sea  breeze  fronts  (for 
example,  see  Purdem  [74]  and  Scofield  and  V.'eiss  [76]).  The  convergence  of 
the  daytime  winds  blowing  from  the  sea  and  the  Gulf  toward  the  Florida  penin- 
sula was  one  of  the  factors  sungested  by  Edinger  [67]  which  accounted  for 
Florida's  summer  thunderstorms. 


? . 1 3 . Speerl_-o  IJ  _ Movr^rent 

An  imi>ortant  feature  for  estiiiiation  of  the  severity  of  windstorms  is 
the  spced-of-cel 1 movement.  Displacements  can  be  measured  from  pairs  of  en- 
hanced infrared  images  using  systems  such  as  the  ISI  (Interpretation  Systems 
Incorporated,  "odd  dSO)  available  at  HESS  (National  Environmental  Satellite 
Service)  by  locating  tlie  intersection  of  two  cursors  over  the  initial  and 
final  positions  of  the  ceil.  Displacement  of  arc  clouds  which  represented 
the  leading  edge  of  a gust  front  was  compared  by  Gurka  [76]  with  low-level 
wind  speed.  The  results  of  the  analysis  revealed  that  "rapidly  moving  arcs 
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arc  generally  associated  with  strong  low-level  winds". 


2.14.  Me  rq j_n  g_C  ej  1 s 

Investigations  by  Simpson  and  Woodley  [71]  indicate  that  merging  cumu- 
lonimbus cells  "often  produce  more  than  an  order  of  magnitude  more  rain  than 
isolated  clouds  on  the  same  day,  probably  owing  to  dynamic  invigoration  of 
the  merged  cloud  circulations".  The  authors  suggested  that  "merger  and 
organization  are  probably  the  first  necessary  steps  in  the  formation  of 
squall  lines,  tropical  storm  rainbands,  and  the  giant  cumulonimbus  systems 
that  fuel  the  large-scale  equatorial  air  motions".  Procedures  for  increasing 
the  estimate  of  convective  rainfall  based  on  occurrence  of  merging  thunder- 
storm cells  were  incorporated  into  Scofield's  model  for  convective  rainfall 
estimation  (Fig.l).  The  merger  of  thunderstorm  cells  can  be  clearly  seen  in 
sequences  of  enhanced  infrared  images.  Merger  can  also  be  detected  in  se- 
quences of  high-resolution  visible  imagery. 

2.15.  Splitting  Cells 

Splitting  cells  and  right  (or  left) deviating  stornis  are  two  severe  storm 
features  which  usually  occur  together.  A cumulonimbus  cell  is  split  into 
two  separate  sections  when  one  of  its  members  is  rotating.  The  cyclonically 
rotating  inoiaber  moves  to  the  right  (right  deviating  storm)  and  tends  to  be 
tornadic.  The  rigfit  moving,  tornado  producing  Davis  storm  studied  by  Burgess 
and  Brown  [73]  resulted  from  the  splitting  on  April  19,  1972,  of  an  iso- 
lated thunderstorm  over  Western  Oklahoma.  After  the  split,  the  left  moving 
section  traveled  2?  degrees  to  the  left  of  the  mean  tropospheric  wind  while 
the  Davis  storm  (the  right  moving  section)  traveled  25  degrees  to  the  right 
of  the  mean  wind.  Figures  5-F-6  (i  through  1)  of  Anderson  [74]  illustrate 
the  correspondence  between  splitting  cells  and  tornado  threat  areas  in  both 
radar  data  and  visible  satellite  dat.«.  Splitting  cells  can  also  be  observed  in 

sequences  of  enhanced  infrared  imaciciy. 
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3 . ^or  Automatic  Detection  ,ind  Trtickinn  of  Severe  Storm  Features 

The  algorithms  presented  in  this  section  address  tne  four  major  problem 
areas  belov;; 

1)  cloud  classification 

2)  cloud  segmentation 

3)  cloud  geometry 

4)  cloud  tracking 

A reviev;  of  algorithms  applicable  to  the  above  topics  and  suggestions  for  future 
algorithmic  development  are  included  in  this  section.  Cloud  classification  is 
the  first  step  in  identification  of  severe  storm  features.  In  order  to  deter- 
mine wind  fields,  clouds  must  be  classified  into  height  categories.  A descrip- 
tion of  cloud  type  is  implicit  in  the  nomenclature  itself  of  some  of  the  severe 
storm  features.  Low-level  cumulus  cloud  line  feeder  bands  consist  of  clouds  of 
cumulus  cloud  typo.  Merging  cells,  splitting  cells,  etc.,  are  features  of 
cumulonimbus  clouds.  Positive  vorticity  centers  are  best  defined  as  middle- 
level  alto  clouds  (Clark  [7^).  When  using  GOES  imagery  in  severe  weather  fore- 
casting, observing  cloud  types  and  their  change  with  time  is  used  to  infer 
changes  in  thermodynamic  instability  (Purdom  [76b]).  If  cloud  classification 
systems  classify  on  an  area  basis  instead  of  a point-by-point  basis,  cloud- 
segmentation  tccliniqucs  must  be  applied  to  obtain  an  outline  of  cloud  features. 
The  cloud  segmentation  of  an  area  of  satellite  imagery  can  be  partial  (looking 
only  for  cloud  features  of ’a  desired  cloud  type  and  shape)  or  total.  If  the 
cloud  segmentation  algorithms  result  in  outlining  and  labeling  cloud  objects 
by  cloud  type  only,  shape  properties  (such  as  comma-tail,  arc-shaped,  line- 
shaped) must  bo  determined  by  algorithms  designed  to  characterize  the  geometrical 
properties  of  the  segmented  cloud  objects.  To  detect  on-going  changes  in  shape 
and/or  motion,  cloud  objects  must  be  tracked  from  frame  to  frame. 
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3.1.  Cloud  Cla^'.r.i f icotion 


A brief  survey  of  cloud  classification  studios  is  followed  in  Sections 
3.1.1.  and  3.1.2.  by  a more  coiiiprohensi ve  discussion  of  major  cloud  classi- 
fication studies  arranncd  witliin  each  section  in  chronological  order.  Methods 
to  discriminate  clouds  from  clo'r  areas  are  presented  in  Section  3.1.1.  The 
developmotit  of  autom.'. t.ic  cloud-type  classification  systems  is  reviewed  in  Sec- 
tion 3.1.2. 

The  problem  of  distinguishing  clouds  from  clear  areas  was  apnroachod  either 
from  the  standpoint  of  selection  of  approp'riate  thresholds  or  weighting  func- 
tions. Arkinn  [6Aj,  f;atz  [6d],  and  Stamm  and  Vonder  liaar  [70]  investigated 
methods  for  selection  of  a single  threshold  from  visible  measurements.  These 
methods  included  manual  selection  of  thresholds  from  series  of  thresholded 
images  and/or  comoca risers  line  plots  of  relative  radiance  vs.  samole  number 
and  autonsatic  selection  of  thresholds  from  gradient  images  and  histograms  of 
visible  data  within  a grid  area.,  Shenk  and  Salomonson  [71]  explored  the  ef- 
fect of  varying  the  ratio  R of  areal-cloud-sizc-to-arcal  resolution  element 
on  estin’.atos  of  actual  cloud  cover.  They  concluded  that  the  value  of  R must 
be  at  least  IQO  if  cloud-cover  estimates  obtained  by  single-threshold  methods 
were  to  bo  valid,  f’otit  Shen!;  and  Salomonson  [71]  and  Miller  and  Feddes  [71] 
designed  ■'v.iltipio  tin'cshold  i”cthods  using  weighting  functions  to  estimate  per- 
centage of  cloud  cover. 

Most  of  the  cloud  classification  studies  were  conducted  using  cloud-type 
information  over  oceans.  Over  the  oceans,  the  clear  area  is,  in  general, 
homogeneous:  radiance  spread  is  less  than  over  land  and  difference  in  radiance 
I levels  between  ocean  and  clouds  is  greater  than  that  hetween  land  and  clouds. 
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Threshold  selection  techniques  to  determine  see-surface  temperature  were 
developed  by  Smith  ct  al.  [70],  Smith  and  Rao  [72],  and  Loose  et  al.  [71].  In 
tropical  oceanic  areas,  cither  the  sea-surface  temperature  (or  a linear  func- 
tion of  the  sea-surface  temperature)  could  be  used  to  discriminate  clouds  from 
background.  The  linear  function  used  by  KoffT'i' et  al,  [73]  to  separate  clouds 
from  clear  area  was  T - S'^C,  whore  T was  the  sea-surface  temperature  for  the 
given  grid  area. 

Cloud-type  classification  systems  using  only  visible  data  characterized 
cloud  types  by  form,  pattern,  texture,  and  dimensions  or  layers  of  the  patterns 
and  forms.  The  guides  to  interpretation  of  TIROS  satellite  pictures  prepared 
by  Conover  [C2]  motivated  several  of  the  earlier  studies.  Rosenfeld  [65]  de- 
termiined  optimal  window  sizes  for  discrimination  of  cloud  patterns  into  cate- 
gories such  as  "brokenness",  elongation,  fi’brosity,  and  convexity.  Leese  [64] 
used  statistical  methods  to  determine  synoptic-scale  features  for  classifying 
low-level  cumuliform  cloud  patterns  which  occur  over  the  oceans  by  cloud  amount, 
size  of  the  cumuli for:ii  cells,  and  size  of  the  cloud  bands.  Cell  size  was 
found  to  bo  related  to  the  location  of  the  observation  and  the  wind  velocity 
field.  Cloud  amount  was  related  to  tlie  location  of  the  observation  with  ref- 
erence to  t!io  center  of  the  anticyclone,  to  tfie  latent  heat  transfer  between 
the  ocean  surface  and  ilie  atmosphei’o,  to  the  wind  fields  at  the  surface  and 
at  tlie  850  ml'  level,  and  to  the  air-dewpoint  temperature  difference  at  the 
surface.  Surface  wind  speed  wms  also  an  important  factor  in  discriminating  * 
among  groups  of  cloud  bands.  Katz  derived  motion-invariant  features  based  on 
the  autocorrelation  function  to  measure  the  direction  and  the  degree  of  align- 
ment of  cloud  streets.  Darling  and  Josepii  [68]  tested  seven  different  pattern 
recognition  techniques  for  discrimination  of  noncuirulus  clouds  from  cumulus 
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clouds  and  solid  cell  cumulus  clouds  from  |iolyc]orial  cell  cumulus  clouds. 

The  screenituj  isultiple  rofjrcssion  technique  with  four  or  fewer  discriminators 
performed  as  \/ell  on  i fid'-pendcnt  samples  as  other  techniques  which  used  300  to 
TOO  properties.  Cootii  [73]  used  multiple  discriminant  analysis  techniques  and 
perceptron  training  methods  to  isolate  sinql e- layered,  low-level  cloud  regions. 
In  order  to  resolve  ambiguities  between  cloud  classes,  tfie  author  reconimendcd 
tne  addition  of  infrared  features  to  the  pattern  vector. 

'with  the  availability  of  nul  t i spectra  1 satellite  data,  automatic  systems 
for  d iscrim i ra tion  of  cloud  types  (as  opposed  to  cloud  forms  or  shapes)  were 
developed.  Koffler  eL  al.  [73]  used  ("ational  Meteorological  Center)  tem- 
pcratin'o  values  \o  classify  clouds  by  height.  Greaves  and  Chang  [70]  deter- 
mined clo'.;''-iypo  sigi:atures  from  dcsi';ni  samnles  of  l,T':'dS-3  cloud  data  over  the 
continental  "i;i  ted  Slaios  in  whicii  tlic're  was  relatively  uniform  cloud  cover  and 
no  raoid  cie  -'inn  or  buildup.  The  point-by-point  classification  mot'nod  of 
Greaves  an-'  Crwang  [70]  suocossfu  1 ly  identified  all  cloud  types  except  for  cumu- 
lus and  cu' sulonir  bus.  T!ic  lad;  of  success  for  cumulus-type  clouds  was  attrib- 
uted to  the  nonhonioqeneous  nature  of  the  cloud  field.  A combination  of  thres- 
holds dorivoT  fron;  a radiative  tramsfer  model,  climatological  knowledge  of  the 
uppe>'  boundc'-ies  of  cloud-tyi'e  surfaces  for  a given  area,  and  expected  roflcct- 
cance  valu>:-  f.'r  cloud  typos  was  used  by  Shciik  and  Holub  [73]  ond  Shenk  et  al. 

[ 7f'j  to  ('utul;lish  four-channel  cloud-type  sionatures  for  cloud  types  over  tropi- 
cal oceans.  Cocth  [73]  combined  statistical  pattern  recognition  techniques 
with  icature  cx*-raction  metliods  to  classify  higli-resolution  , dual-channel 
[;0AA-1  cloud  samples  over  tropical  oceanic  regions.  A comparative  study  of 


selected  feature  sets  and  s'atistical  pattern  recognition  techniques  to  classi- 
fy tropical  oceanic-cli  tyni^  - was  conducted  hy  Parikh  [77a].  Classification 


32 


accuracy  deteriorated  depending  on  the  complexity  of  the  cloud-type  decision 
problem.  Optimal  thresholds,  either  for  point-by-point  cloud-type  signatures 
or  for  decision  regions  in  feature  vector  space,  did  not  remain  constant  v/ith 
change  in  location  or  season  of  observation.  A method  of  cloud  classification 
by  cloud-sample  segmentation  and  comparison  of  segment  features  was  developed 
by  Parikh  [77b].  The  method  was  independent  of  design  samples  or  prior  thres- 
hold selection  techniques.  All  of  the  above  mul tispectral  techniques  except 
for  that  of  Greaves  and  Chang  [70]  have  been  tested  primarily  or  exclusively 
on  cloud  data  over  oceanic  regions.  Further  study  is  needed  to  determine  what 
modifications  must  be  made  cither  in  cloud-type  signature  classification  sys- 
tems or  feature-extraction  classification  systems  to  discriminate  cloud  types 
over  land. 
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3.I.]  . for  f)i 1_ i r.'"! u i r.h i !Ki  Clciiid",  Fron 

c:jc;'ir  /.rcj  i:; 


v.’itli  bricjhtnosr,  valuer;  above  c were  dcoignn  tc.d  as  clcuds  and 
point::  v.'ith  brightness  values  equal  to  or  less  than  c v.-are 
designated  as  bachground.  Except  for  snox'-  cr  ico-coverod 
regioris,  regions  cf  thin  cirrus,  or  rogienn  of  sr.all  and 
scattered  convective  cloud  elcrrcnts  belov,'  sensor  resolution 
sine,  a brichtress  threshold  should  be  a good  criterion  for 
the  detection  of  the  presence  cf  clouds  acainot  a carl;  bac?.- 


crour.;:  suen  as  tne  ccea: 


face.  hccordinn  tc  Starr,  an: 


, a histograr,  of  the  frequency 


cJt,  a 


aiven  radiance  level  occurs  vithin  a grid  area  will  contain 
a peal:  cf  rouchly  Gaussian  shape  representing  a clear  area 
provider!  the  clear  area  is  hor.cgcnccus  (ocean  surface  rathe: 
than  land  surfa:-c).  Hcwc;vcr,  even  fer  ocean  surfaces  it  is 
difficult  to  detcrrinc-  tlic  f.'oint  tc  tlio  riec.L:  of  the 
Gn,us:;i  an  C-loar-arca  iicn!:  at  v.'hich  clc'vv.l  interference  is  bo- 
c oin i n ■ .i  npor t,  a 1 1 1 . 


Limitations  of  son.sor  -systems  us  well  ns  variable 
briglitnes.s  re.sponses  for  a cloud  of  given  reflectivity  must 
be  considfrred  for  both  mum:  1 and  nutomatic  estimntes  of 


cloud  cover.  T)iin  cirruu  ir.  tranr.p:ir(;nt  or  scmiLrnnsparcnL 
and  Uic  reflectance  characteriuticn  of  the  underlying  back- 
ground detcniiine  the  recorded  brightnen.o  intensity  level, 
fmall  cumulu:;  c]ov;dn  below  rf? solution  cannot  be  oeen 

on  t)ic  image  );ut  tlic-ir  bri.ghtnene  affects  tlic  sensor  reading. 
Variabilities  of  tlie  signal  witliin  a camera  system,  changes 
of  response  from  one  system  to  another,  weakness  in  the  cross 
track  normalization  technique,  changes  in  camera  calibration 
and  degradation  in  camera  response  substantially  alter  the 
brightness  response  of  individual  cloud  elements  from  day  to 
day . 

In  19&4 , Arking  [Gi]  discussed  the  problem  of  cloud/back- 
ground threshold  selection  in  connection  with  automatic 
analysis  of  Tiros  III  television  pictures  to  determine  lati- 
tudinal distribution  of  cloud  cover,  /irking 's  approach  v.'qs 
to  employ  the  subjective  judgment  of  a hum.an  observer  to 
compare  a scries  of  two-level  throsholdcd  images  to  the 
oriqin,.!  .••.atelbite  image  to  determine  which  of  the  throsholdcd 
iinaqcs  no;-t  faithfully  reproduced  tlie  c]oud  boundaries  of  the 
original  inriagc.  For  each  orbital  sequence  a new  threshold 
was  chosen.  It  was  found,  however,  that  the  same  threshold 
could  bf'  used  for  each  satellite  picture  v.’ithin  a given  orbit. 

Ka  tr,  constructed  a gradient  image  to  automatically 
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determine  a brightness  throshoid.  The  gradient  image  is  de- 
fined as  the  image  whose  value  at  every  point  is  given  by 

II  G II  = (dl/dy)^ 

v.’hcre  1 is  the  brightness  value  of  the  given  point.  ! I G | | 
“"an  be  approximated  by  ] ] G | | w.here 

II  G II  = max{  I 3I/8x 1 , j ?I/cy I } + j min { ! ? I/oX | , 1 S I/3y | } . 

In  order  to  facilitate  handling  the  pictorial  data  in  packed 
form.,  the  quantity  actuf^lly  used  for  com.putation  of  the  gra- 
dient im'.age  was  2 ( |g|  | . The  cumulative  frequency  percentile 
was  calculated  for  each  gradient  level  in  the  gradient  image. 
The  smallest  gradient  level  g for  which  the  cumulative  fre- 
quency percentile  v;as  greater  than  or  equal  to  a given  value 
p was  then  selected.  Expcrim.ental  results  vrerc  calculated 
for  p = 75%,  p = 87.5%,  and  p = 93.7%  v.’ith  the  value  of 
p = 87.5%  the  personal  choice  of  the  author.  The  average 
intensity  of  those  points  at  which  the  gradient  was  at  least 
g was  used  as  the  threshold  for  cloud/background  separation. 

Shenk  and  Galomonson  [7]]  suggested  using  two  bright- 
ness thresholds  to  divide  the  satellite  image  into  three 
data  sets  --  data  v.’here  the  resolution  element  is  clear,  data 
wherr'  the  resolution  element  is  partially  cloud  filled,  and 
data  where  the  resolution  clement  is  completely  cloud  filled. 
Let  be  the  number  of  points  in  the  first  data  set,  N2  the 
number  of  points  in  the  second  data  set,  and  the  number  of 
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points  in  the  third  data  set.  Then  the  percentage  of  cloud 
cover  was  calculated  as 

N„+2N^ 

N = 2 3 

2(Nj^+N2+N^) 

The  association  of  a v;cight  of  1 with  elements  in  the  second 
data  set  and  a weight  of  2 with  elements  in  the  third  data 
set  implies  that  a cloud  cover  of  50%  is  assumed  for  every 
partially  cloud-filled  resolution  element.  This  two-threshold 
method  is  contrasted  with  the  single- threshold  method  which 
divides  a satellite  image  into  points  w’here  the  resolution 
element  is  clear  and  points  w^here  the  resolution  element  is 
partially  or  completely  cloud  filled.  For  a typical  situation 
in  the  tropics  in  which  the  average  cumuliform  cloud  size  is 
4 km  and  for  a typical  sensor  resolution  of  4 km,  the  single- 
threshold mctliod  would  estimate  a true  cloud  cover  of  20%  as 
70%  cloud  cover.  The  authors  concluded  from  studies  on  simu- 
lated data  and  high-spatial  resolution  (300  feet)  Apollo  6 
photographs  that  the  value  of  the  ratio  R of  areal  cloud  size 
to  areal- resolution  element  must  be  at  least  100  if  cloud 
cover  estimated  from  single- threshold  methods  and  actual 
cloud  cover  are  to  agree  within  approximately  10%.  However, 
if  cloud  cover  is  estimated  by  the  two-threshold  method,  good 
cloud  cover  estimates  result  for  values  R ^ 10.  If  the  value 
of  R is  known,  nomograms  presented  by  the  authors  can  be  used 
to  obtain  true  cloud  cover  from  estimated  cloud  cover. 

The  question  of  automatic  selection  of  the  one  or  two 
brightness  thresholds  did  not  arise  in  the  study  of  Shenk  and 
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Salomonson  [71].  The  single-threshold  method  applied  to 
simulated  and  high-resolution  cloud  data  involved  counting 
manually  for  each  of  seven  resolution  sizes  the  number  of 
resolution  elements  which  contained  no  cloud  elements.  Simi- 
larly, for  the  tv.’o-tlireshold  method,  a manual  count  was  made 
of  those  resolution  elements  which  were  partially  cloud 
filled  and  of  those  v/hich  v;ore  completely  cloud  filled. 

Automation  of  a multiple  threshold,  weighted  histogram 
model  for  cloud-cover  estimation  was  described  by  Miller  [71] 
and  Miller  and  Feddes  [71].  Satellite  image  data  compressed 
from  a 0-63  scale  to  a 0-14  scale  was  divided  into  the  five 
equally-spaced  brightness  ranges  0-2,  3-5,  6-8,  9-11,  and  12- 
14.  With  each  brightness  range  was  associated  an  empirically- 
derived  seasonal  weight.  For  October-May,  the  weights 

, i = 1,...,5,  for  the  five  brightness  ranges  were  0,  2,  7, 

8,  8,  respectively.  For  June-Septem±)er , the  w’eights 

lA  , i = 1,...,5,  for  the  five  brightness  ranges  were  0,  2.5, 

7.5,  8,  8,  respectively.  Higher  weighting  of  classes  2 and 

3 for  tlie  summer  months  reflected  the  increased  frequency  of 
occurrence  of  small  cumulus  cloudiness  during  the  summer. 

The  percentage  K of  cloud  cover  was  calculated  as 

f W.F. 
i=l  ^ ^ 

N = ( -i-± ) + k 

64 


W’here 

IvA  = weighting  factor  for  class  i 

= frequency  of  occurrence  of  the  ith  class 

k = 0.51,  a computer  t.uncation  constant. 
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The  values,  although  in  general  agreement  with  concurrent 
surface  observations,  showed  a consistent  tendency  to  under- 
estimate surface  reports  of  cloud  cover.  Miller  and  Feddes 
t ascribed  the  differences  primarily  to  differences  be- 
tween tlie  field  of  view’  of  the  surface  c-bserver  and  the  angle 
of  view  of  the  satellite  cam.cra  (which  resulted  in  an  over- 
estim.ate  of  cloud  cover  from  surface  observations)  and  to  the 
lo’wer  response  threshold  (i.e.,  lack  of  sensitivity  to  thin 
cirrus  or  to  si.iall,  scattered  convective  elcm.cnts)  and  re- 
solving pov.'cr  of  the  satellite  sensor.  The  authors  concluded 
that  the  automated  values  for  were  at  least  as  good  as 
those  which  could  be  derived  by  eye  from,  satellite  pictures. 

A threshold  T can  be  applied  to  infrared  m.easurcmcnts  of 
cloud  data  to  discriminate  clouds  from  background.  Over  the 
occiins,  an  ap;.>ropr iatc  value  for  f for  each  grid  area  of  the 
image  would  be  cither  the  sea  surface  tcm.perat’ure  T within 
that  grid  area  or  T-k  where  k is  some  sm.all  positive  constant. 

Kofflor  ct  al.  [73]  subtracted  a constant  of  5°C  from  the  sea 
surface  tem.pcrature  to  determine  the  grid  area  thresholds  T. 

/ill  grid  points  with  infrared  tcm.peraturc  greater  than  or 
equal  to  T were  classified  as  background  points. 

A statistical  histogram  method  to  determine  the  sea  sur- 
face temperature  values  T from  niglittimc  3.8  pm  infrared  data 
is  given  in  Smith  ct  al.  [70].  Tlic  following  procedure  was 
used  to  analyze  the  temperature  histogram.  If  a clear  dis- 
tribution (Gaussian)  could  be  detected  such  that  the  maximum  fre- 
quency of  a point  of  the  clear  clist  rJ  !)’ution  war.  both  at  a tcmpcratv;rc- 

39 


. . .. 


greater  than  freezing  {273°k)  and  greater  than  10  percent  of 
the  total  number  of  observations,  then  the  histogram  was  con- 
sidered to  be  sufficiently  cloud  free  to  be  acceptable  for 


f 


i 


u 


the  second  stage  of  analysis.  Find  the  point  ^g^ope  ^ 

high  temperature  wing  of  the  clear  distribution  where  the 
change  of  frequency  with  temperature  is  a maximum, i . e . , the  I 

point  of  maximum  slope.  Define  the  highest  observed 

temperature  value  with  a frequency  greater  than  1 percent  of 
the  total  number  of  observations.  Define  T as 


where 


T 


M 


slope 


o 


0 H known  standard  error  of  measurement 
For  this  study  (Nimbus  2),  a = 1.5'’k.  T will  usually  be 
greater  than  the  modal  peak  temperature  of  the  clear  distri- 
bution when  clouds  exist.  If  the  difference  between  M,  . , 

high 

and  T is  greater  than  3a  or  if  the  change  of  frequency  v?ith 
temperature  at  ^'^g^ooe  less  than  3 percent  per  degree  K, 
then  the  high  temperature  wing  of  the  clear  distribution  is 
influenced  by  cloud- contaminated  measurements  and  sea- 
surface  temperature  cannot  be  specified.  Otherwise,  specify 
the  sea-surface  temperature  as  T. 

Sea-surface  temperatures  for  areas  of  severe  cloud  con- 
tamination were  inferred  through  space  and  time  interpola- 
tion procedures.  In  a few  cases  in  which  the  histogram 
technique  failed  to  differentiate  between  cloud-free  and 
cloud-contaminated  temperatures,  errors  were  attributed  to 

the  existence  of  low-level,  uniform,  overcast  conditions. 
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Tlic  authors  suggested  using  either  spatial  gradient  consis- 
tency chocks  and/or  time  compositing  techniques  to  filter  out 
erroneous  sea  temperatures.  Relative  discrepancies  between 
the  histogram  sea-surface  temperatures  and  ship-reported  sea- 
surface  temperatures  were  generally  less  than  1°K. 

Smith  and  Rao  [ 72J  developed  an  algorithm  to  determine 
both  daytime  and  nighttime  sea-surface  temperatures  from 
spatially  contiguous  sets  of  simultaneous  3.7ym  and  Hum 
infrared  measurements.  The  contiguous  resolution  elements 
were  assur.ed  to  have  the  same  average  cloud  tem.perature  and 
sea-surface  temperature,  i.c.,  the  same  type  of  cloudiness 
co.nditions.  Fcr  resolution  elem.ents  1 and  2 and  spectral 
windov.’s  3.7  ym  and  11  ym,  the  ratio  of  the  cloud  am.ounts  of 
the  two  resolution  elements  v:ere  equated  for  each  spectral 
window,  i.e., 

^3.7  ~ ^3.7 
^3.7  ■ ^3.7 

whore 

I = miCasured  radiance  for  resolution  element  r 

for  spectral  wavenumber  S 

B (T_)  = rianck  radiance  for  spectral  v.’avenumber  S 

O ^ 

The  Planck  radiance  B_(T)  is  a function  of  temperature  which 
is  given  by 

Bs(T)  = (BXP  (C2S/T)  - 1) 

where 


hi  - hi'h> 


hi  - hi'h> 
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r 


= 1.19061  X 10 


-1 


erg-cm  -sec 

C2  = 1.43868  cni-deg 

If  the  relation  R(T)  is  used  to  denote  the  result  of  equating 
the  ratio  of  cloud  amounts  for  the  two  spectral  windows,  i.e., 


R(T) 


can  solve  for  Tg  (sea -surface  temperature)  by  settinq 
R(T)  = 0.  The  first  order  Taylor  approximation  for  T (where 
the  superscript  k denotes  a quantity  evaulated  for  the  kth 
iteration)  is  given  by 

^k+1  ^ ^k  _ jpk/(  ?R  jkj 


where 

?R  k _ 2 _ 1 ^ _ 2 

V 3T  ^-^3,7  3.7^  ^ 9T  ’ ^"11 


1 ^^3  7 

I ] ( 
tv' 


£T 


vrith  an  intial  guess  temperature  of  310°k,  a solution  was 
usually  obtained  in  less  than  five  iterations.  This  multi- 
spectral  approach  was  one  met}iod  used  by  Smith  to  resolve 
the  problemiS  of  interpretation  of  daytime  3.7  tm  data  in 
whic}i  cloud  contaminated  temperatures  can  be  equal  to  or 
larger  than  clear  sky  brightness  temperatures. 

Another  approach  suggested  by  Smith  for  determining  sea- 
surface  temperatures  using  single- channel  data  was  described 
in  a paper  by  Loose  et  al.  [71].  This  method  assumed  a 
Gaussian  distribution  of  the  system  noise,  but  unlike  the 
method  of  Smith  et  al.  [70]  previously  reviewed,  the  value  0 
of  the  standard  error  of  measurement  (noise)  was  not  assumed 
to  be  known.  Given  any  three  temperature  classes  of  the 
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clear  distribution,  the  sea-surface  temperature  Tg  which  is 
the  mean  value  of  the  distribution  from  clear  atmosphere 
sajnples  was  computed  as 


T = 
S 


(T3^  - T3^)ln  F3/F2  + - T2^)ln 

2(Tj  - T^)  In  F^/F^  - 2(T^^  - T3)ln  F^^/F^ 


where 

, i = 1,  2,  3 is  the  ith  temperature  class 
F^ , i = 1,  2,  3 is  the  frequency  of  occurrence  of 
tem.perature  class  T^. 

A discribution  of  values  of  Tg  was  obtained  by  using  a number 
of  ccm±)ination5  of  temperature  classes  taken  three  at  a 
time.  The  mean  value  of  this  second  distribution  was  used  as 
an  estimator  for  the  sea-surface  temperature  Tg.  Six  cloud 
and  noise  contamination  tests  had  to  be  passed  by  the  second 
distribution  of  values  for  Tg  for  the  sample  to  be  retained 
for  sea-surface  temperature  calculations.  Since  cloud  layers 
can  produce  histograms  of  the  saime  shape  as  those,  obtained 
from  a clear  atmosphere,  Leesc  suggested  that  instead  of  just 
restricting  Tg  to  be  above  273°k,  a suitable  test  should  be 
developed  on  the  basis  of  first-guess  fields.  A first-guess 
field  for  the  sea-surface  tem.perature  value  could,  for  ex- 
ample, consist  of  the  previous  day's  temperature  combined 
with  appropriate  limits  for  acceptable  twenty-four  hour 
changes  (obtained  from  gradient  calculations).  Results  in- 
dicated that  errors  of  less  than  1°C  should  be  possible  for 
areas  where  the  temperature  gradient  is  less  than  2°C  per 
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100  km.  For  reyions  wheic  tho  tomi a Lur e gradient  was  2 to 
4^0  per  100  km,  visual  data  should  be  included  to  reduce 
cloud  contamination  errors.  Improvcmonts  may  be  needed  in 
earth  location  procedures  if  an  enor  of  less  than  1°C  is  to 
be  maintained  for  regions  in  wliich  the  temj:>eraturo  gradient 
is  larger  than  4°C  per  100  km. 

Stamm  and  Vender  Haar  [70  J sugejested  <is  a result 
of  threshold  studies  in  the  visible  channel  that  some  com- 
bination of  threshold  criteria  in  several  spectral  channels 
would  be  able  to  detect  cloud  contamination  with  better  pre-  ; 

cision  than  si ngle-cha.nnel  thresholds.  They  investigated  two 
threshold  methods  for  the  visible  cliannel.  The  first  method 
compared  a line  plot  of  relative  radiance  vs.  sample  number 
to  a corresponding  line  on  the  given  picture  to  manually  de- 
tex'mine  a threshold  between  cloud  and  background  regions. 

The  second  method  analyzed  the  histograms  of  radiance  levels 

v.'ithin  a given  grid.  For  grid  areas  over  oceans  with  a i 

reasonable  amount  of  clear  area,  a large  peak  of  roughly 

Gaussian  shape  represented  the  clear  area.  By  forming  a 

Gaussian  curve  from  the  data  on  the  left-hand  side  of  the 

peak  (the  side  free  from  cloud  interference) , it  was  possible 

I 

to  predict  the  clear  area  peak.  The  point  at  which  the 
Gaussian  curve  starts  to  rise  again,  or  a few  radiance  levels 
to  the  cloudy  (right)  side  of  this  clear  area  peak,  was 

chosen  as  the  threshold.  \7hcn  in  doubt,  the  threshold  was  i 

chosen  to  overestimate  the  cloudy  regions. 
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3.1.2.  Automatic  Cloud-Typo  Classification  Systems 

Various  automatic  cloud  classification  systems  using 
features  derived  from  visual  satellite  data,  infrared  satel- 
lite data,  and  multispectral  satellite  data  have  been  devel- 
oped from  1964  to  the  present  date.  In  this  section,  a brief 
survey  of  techniques  and  results  of  single-channel  cloud  clas- 
sification systems  is  follov.'od  by  a more  comprehensive  review 
of  mul  ti  sijoctral  systems. 

Conover  [c,2  ] published  in  1962  a guide  to  cloud  inter- 
pretation from  satellite  pictures.  Features  selected  to  dis- 
criminate cloud  typos  were  form  (round,  curved,  or  elliptical), 
pattern  (banded  or  randomly  spaced)  , te.xture  (smooth  versus 
fibrous),  brightness,  structure,  and  dimensions  of  the 
patterns  and  form.s.  Each  cloud  type  was  associated  v;ith  one 
of  four  brightness  classes  --  dark  gray,  gray,  white,  and 
very  white.  The  author  suggested  use  of  automatic  miethods 
for  cloud  reflectivity  contouring  as  an  aid  to  photointerpre- 
tation of  satellite  images.  The  idea  of  characterizing 
cloud  types  by  ranges  of  measurements  in  one  spectral  window 
was  a forerunner  to  completely  automated  cloud-type  classifi- 
cation schemes  such  as  that  of  Shenk  and  Holub  f72]  in  which 
ranges  of  measurements  in  four  spectral  windov.-'s  specified  a 
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unique  cloud  type.  j 

In  19G4,  Leese  [G4]  applied  rnultiple-discrirainant 
analysis  tec’nniqucs  to  tlic  problein  of  classifying  low-level  ’ 

curr.ulif orr;.  cloud  patterns  by  cloud  amount,  size  of  the  cum.u- 
lifcr.’u  cells,  and  size  of  the  cloud  bands.  Tv, 'cnty- three 
synoptic  features  consisting  of  variables  such  as  latitude, 
longitude,  mean  solar  tine,  surf ncc-gecstrophic  wind  speed, 

, 

and  latent  heat  flux  v;ore  used  to  categorize  the  cloud 
patterns  into  cloud  cover  groups,  cloud  size  groups,  and 

cloud  banding  groups.  Cloud  amount  was  found  to  be  most  • | 

1 

closely  related  to  the  location  of  the  observation  with  re-  i 

ference  to  the  center  of  t.he  anticyclone.  Size  of  the 
cum.ulus  cells  v.'as  most  closely  related  to  the  structure  of 
the  wind  field  between  the  850-mb  level  and  the  surface. 

Discrimination  results  betv.’eon  cloud  banding  groups  were  not  . 

as  successful  as  discrimination  results  betv.’een  cloud  am.ount  ■ 

groups  and  cloud  sizi  groups.  More  objective  techniques  for  ’ 

the  original  description  of  the  cloud  patterns  (such  as  in- 
creased number  of  cloud  amount  groups,  information  on  cloud 
patterns  on  different  size  scales,  and  a coordinate  system 
to  describe  the  location  of  the  observation  in  the  anti- 
cyclone) combiined  with  improvements  in  the  sensor  system  of 
the  satellite  v.’cre  suggested  as  possible  factors  which  would 
improve  classification  results.  In  particular,  two- 
dimicnsional  spectral  analysis  teclmiques  such  as  numerical 
filtering  should  be  used  to  pjrovide  information  about  cloud 
patterns  on  at  least  three  scrales  --  tlie  synoptic  scale 


(more  than  50  mi) , the  mcsoscale  (between  10  and  50  mi) , and 
the  submesoncale  (less  than  10  mi) . 

Rosenfeld  [65]  discussed  a scanning  v;indov/  approach  to 
automatic  clas.sif ication  of  cloud  patterns  into  categories 
representing  r.parseness  or  "brokenness " , fibrosity,  elonga- 
tion, regularity,  straightness,  and  convexity  properties. 
Quantitative  measurement  scales  for  these  parameters  v;ere 
derived  from  psychological  studies  of  human  judgments  of 
these  parameters  in  a series  of  cloud  photographs.  In  order 
to  discrim.inate  between  solid  black  (background)  , solid 
white  (cloud),  and  broken  cloud  regions,  the  number  of  back- 
ground elements  in  the  scanning  square  window  of  an  im,age 
was  counted  for  each  position  of  the  window.  If  the  number 

of  backerreund  elements  fell  below  a threshold  T , the  scuare 

w 

was  cailcd  solid  white.  If  the  number  of  background  elem.ents 
was  greater  than  a threshold  T,.,  where  a T , the  scuare  was 
categorised  as  solid  black.  Otherwise,  it  was  considered  as 
broken  cloud.  Boundaries  between  regions  of  different  types 
.vere  m.arkcd.  Suggested  window  sizes  for  discrimination  of 
aloud  patterns  with  respect  to  the  above  parameters  were 
5x5  resolution  elem.ents  for  "brokenness"  and  15x15  resolution 
elements  for  shape  and  size  properties. 

ICatz  [64  ] derived  m.otion-invariant  features  to  describe 

the  degree  and  direction  of  alignment  of  striated  clouds  and 

considered  features  which  would  give  inform, ation  about  the 

distribution  of  cloud  sizes  in  a region.  Tiio  ratio  of  the 

axes  of  the  ellipse  of  concentration  of  tlie  central  hump  of 

the  autocorrelation  function  p(li,k)  was  used  to  measure  the 
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degree  of  alignment  of  cloud  streets.  The  direction  of  the 
major  axis  was  taken  to  bo  the  direction  of  the  streets.  The 
ellipse  of  concentration  was  taken  as  the  contour  enclosing 
the  central  hump  of  p(h,k)  for  which  p(h,k)  = 0.2.  The  di- 
rection of  the  major  axis  of  the  ellipse  was  given  by 

2y  11 

tan  20  = — 

^20  ^02 

where  are  central  moments.  The  ratio  r of  the  -axes  is 
r = t where  and  ^2  the  roots  of 

^'^20^02~'^11  ^ ■ ^^*02  ^20^^  + 1 = 0. 

The  autocorrelation  function  p(h,k)  is  the  normalized  (divided 
by  the  individual  variances)  autocovariance  function  C(h,k) 
where 


C (h,k) 


I [I (x,y) -m3 [I (x-h,y-k) -m] dxdy 


with 


I(x,y)  = brightness  level  at  point  (x,y) 
m = mean  value  of  1 


The  feature  a(t)  described  by  Blum  [G23  where  a(t)  is  the 
area  at  time  t of  the  region  enclosed  by  a wavefront  moving 
from  a pattern  at  constant  speed  in  all  directions  was  pro- 
posed to  characterize  the  distribution  of  cloud  sizes.  As 
the  dark  areas  overrun  the  cloud  areas,  the  clouds  begin  to 
disappear  — for  a convex  cloud,  this  would  happen  in  the 
time  it  takes  the  wave  to  traverse  )ialf  the  width  of  the 
cloud . 

Darling  and  Joseph  l(,  l ] applied  five  linear  decision 
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tcchniqucj: , one  nonlinenr  decision  technique  ("Madeline"), 
and  n screening  multiple  regression  technique  to  categorize 
cloud  patterns  into  cumulus  and  noncumulus  classes  and  to 
categorize  curaulus  clouds  into  cumulus  polygonal  cell 
patterns  and  cumulus  solid  cell  patterns.  For  the  linear 
decision  Function  te'chniques  and  for  Madeline,  a group  of  n 
tlireshold  units  bj^(j),  i - l,...,n  vere  randomly  connected 
to  elements  in  the  pattern  space.  The  output  of  the  ith 
logic  unit  b^(j)  for  the  jth  pattern  is  1 if  a function  of 
its  inputs  cMcecds  a threshold  and  0 otherwise.  The  thres- 
hold was  defined  by  a quadratic  switching  surface 


x’’(V 


-1 


-1. 


Ty 


-1  T 

M1-M2 


-1 


M^  + ln  ( 


L\/\U\)=o 


wl'.cre  vectors  M^  and  are  sample  moans  and  and 
variance  matrices  for  classes  1 and  2 respectively  and  X is  the 
pattern  vector.  'Die  form  of  the  decision  function  D(j)  for 
tlie  linear-decision  techniques  v;as 


D (j) 


n 


1 vmb  (j) 
i=l 


6 


w^,i  = l,...,n  are  property  weights 
C is  t!ie  response  unit  threshold 
heights  and  respons'^  function  wore  calculated  according  to 
tin-.'  following  five  algorithms  — forced  learning,  Bayes 
weiqlits,  error  correction,  iterative  design,  and  mean  square 
error.  The  first  throe  methods  are  associated  with  the  a 
porceptron  of  Rosenblatt.  For  Madaline,  a second  layer  of 
threshold-logic  unit?,  was  connected  between  tlie  first  set  and 
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the  response  unit  (decision  function) . 


I 


For  the  screening  multiple  regression  technique,  twenty- 
eight  features  computed  from  a binary  representation  of  the 
scene  were  screened  at  the  0.01  level.  For  this  technique, 
the  decision  function  was  a linear  multiple  regression  equa- 
tion of  the  form. 


whore 


Xo  = + A^X2  +■ 


rn.yCn^+n^)  in  class  1 

~ I 

t-n^/Cn^+n^)  in  class  2 


^0 

A, 

^‘0 

is 

n • 
1 

is 

X • 

is 

■+  A X 
P P 


both  classes 

A^'s  arc  constants  to  be  determined  by  least  squares. 

At  each  step,  the  particular  x^  having  the  highest  correla- 
tion rQ^  v.’ith  Xq  is  added  to  the  equation  if  the  reduction 

in  variance  is  greater  than  0.01.  The  correlation  Zf.- 

Ox  u 1 

is  defined  as 


..  ri"2 


where 

dj^  = difference  between  the  mean  value  of  x^  in  class  1 

and  its  mean  value  in  class  2 

S..  - sum  of  squares  of  x.  over  both  classes. 

11  •“  a. 

The  value  of  A^  is  given  by 
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A . = 

1 


d . n, n„ 
a 1 2 


After  the  selection  of  each  x^,  the  remaining  x^'s  are  ortho- 
gonaiiced  with  respect  to  the  selected  x^^  to  yield  an  equa- 

/s 

tion  for  x^  with  uncorrelated  terms.  The  decision  rule  is ; if 

y-Q  > - (n^^-n^ ) / (nj^+n^ ) ] , assign  the  pattern  to  class  1. 

Otherv.’ise,  assign  the  pattern  to  class  2 if 
^ 1 

Xq  < - 2 [ (n^-n., ) / (rij^+n^  ) ] . For  the  classification  of  non- 
cumulus  versus  cumulus  clouds,  the  features  chosen  were  the 
80- contour  area  (area  bounded  by  the  Y axis  and  0.80  auto- 
correlation contour  for  X > 0) , the  brightness  variance,  and 
the  relative  frequency  of  cloud  size  1 to  25  (i.e.,  from  1 
to  25  resolution  elements). 

Decision  accuracies  for  the  error  correction  procedure, 
iterative  design,  Madalinc,  and  screening  multiple  regression 
teclinicues  were  close  to  901.  The  authors  noted  that  the 
screening  multiple  regression  technique  with  only  three  to 
four  discriminators  performed  as  well  or  better  than  the 
other  tecliniqucs  wliich  used  approximately  400  properties. 

They  recoimcndc'd  enlarging  the  set  of  features  to  include  in- 
foimiation  content  in  the  derivative,  curvature,  gradient, 
and  Laplacian  of  the  brightness  image,  the  coefficients  of 
the  power  spectrum,  and  the  characteristics  of  the  auto- 
correlation function  obtained  by  rotating  a pattern  about  its 
center  of  gravity. 

Booth  [73]  categorized  32x32  visual  arrays  of  ATS-1 
images  into  a 5-group  model  .md  a 3-group  model  using  multiple 
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d i^:criminant  analysis  and  perception  training  methods.  The 
5-grcup  model  consisted  of 

(1)  predominantly  single -layered , lov/-level  8000  ft) 
cloud  regions 

(2)  predominantly  single- layered,  middle-  or  single- 
layered  high-level  cloud  regions 

(3)  predominantly  multi-layered  or  large  convective 
cloud  regions 

(4)  cloud  regions  containing  more  than  one  of  the  above 
groups  vith  one  group  not  dominant 

(5)  none  of  the  above  (i.e.,  overcast  and  clear-sky 
cloud  regions ) . 

The  3-group  model  consisted  of  classes  1,  2,  and  3 above. 
Rcsult'i  of  multiple  discriminant  analysis  using  forty-six 
features  assumed  to  be  normally  distributed  with  the  same 
covariance  matrices  ranged  from  68%  decision  accuracy  for 
group  1 independent  (test)  data  for  the  5-group  model  to 
92%  decision  accuracy  for  group  1 independent  data  for  the 
3-group  model.  Every  niisclassif icd  group  1 sample  was 
labelled  as  group  3 by  the  multiple  discriminant  classifier 
--  thus  pointing  out  the  confusion  between  groups  1 and  3. 
Tests  on  independent  data  v.'orc  performed  exclusively  on 
observations  from  group  1,  since  the  primary  intent  of  this 
classification  system  was  to  isolate  those  regions  which 
would  bo  suitable  input  target  areas  for  the  automated 
operational  system  of  NO/uk  for  computation  of  single -layered , 
low-lcvcl  v.'ind  vectors. 
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Experiments  to  separate  group  1 from  the  combined  groups 
2 and  3 by  pcrccptron  analysis  failed  after  454  iterations. 
The  similarities  bctv.’ecn  groups  1 and  3 obviously  led  to  the 
lack  of  linear  separability  bctv;ecn  observations  from  group  1 
and  observations  from  groups  2 and  3.  rcrccptiron  classifica- 
tion of  787  test  observations  into  cloud  and  no-cloud  groups 
resulted  in  653  (83%)  observations  correctly  classified. 

Suggestions  by  the  author  for  improvem.ent  of  decision 
accuracy  were  to 

(1)  construct  a nev;  5-group  model  consisting  of  low- 
level  single-layered  regions,  middle-high  level 
single-layered  regions,  multi-layered  or  convective 
regions,  overcast  regions,  and  clear-sky  regions, 

(2)  augm.ent  the  pattern  vector  with  features  from  in- 
frared satellite  data  such  as  cloud-top  temperature 
to  help  resolv'c  ambiguities,  for  example,  between 
groups  1 a.nd  3 , and 

(3)  gen'-’ratc  an  optimum  input  pattern  vector  by  per- 
form.ing  scroen.ing  analysis  on  features  derived  from 
the  satellite  data. 

Fatcllite  infrared  radiation  (IR)  tem.pcratures  segmented 
at  each  grid  square  (32x32  IR  spots)  into  four  intervals  de- 
pendent on  synoptic  temperature  fields  obtained  on  a 12-hour 
basis  from  the  National  Meteorological  Center  (NMC)  were  used 
by  Kofflcr  et  al. [73]  to  generate  pictorial  displays  of  cloud- 
type  information  consisting  of  four  gray  sliades.  Each  of  the 
four  gray  levels  represented  a cloud  category  as  shown  in 
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Tabic  2.1  belov,’. 


Gray  Shade 

Temperature  Interval 

Cloud  Type 

Cloud  Top  Jieight 
Interval 

black 

Infrared  Temperature 
a T^-5°C 

no  clouds 

1 

dark  gray 

T^-S^C  > Infrared  Tern- 

low  clouds 

surface  - 700  mb 

perature  > 1’2 

1 

layer 

light  gray 

^ Inf rared  Temper- 

middle 

700  mlD  - 400  mh. 

ature  ^ 

clouds 

layer 

white 

T^  > Infrared  Temper- 

high 

eibove  4 00  mb 

ature 

clouds 

layer 

Tabic  2,1:  Cloud  classification  systen  of  Koffler  et  al.[73].  Threshold 
values  , T^,  T^  represent, respectively,  NMC  temperature 
valuer  per  grid  area  at  the  surface,  700  mb,  and  400  m.b. 


Assumptions  inherent  in  the  use  of  this  classification  scheme 

are : 

(1)  clouds  are  completely  opaque, 

(2)  radiation  em.anates  from  cloud  tops, 

(3)  there  arc  no  rapid  atm.osphcric  changes  in  the  time  in- 
terval bctv:ccn  the  KMC  data  and  the  infrared  temperature 
data , 

(4)  the  earth's  actual  surface  tern, per ature  is  constant  over 
the  area  of  a grid  square  and  equivalent  to  the  NMC- 
supplicd  earth  surface  temperature,  and 
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(5)  errors  caused  by  assignnent  of  an  IR  temperature 
measurement  to  the  lowest  pressure  layer  in  cases 
of  a non-uniquo  correspondence  of  a particular  IR 
temperature  to  a given  pressure  interval  resulting 
from  the  existence  of  an  isothermal  or  near- 
isothermal  atmosphere  or  from  tlic  existence  of  in- 
versions, are  minimal. 

Assumption  number  (1)  leads  to  errors  of  interpretation 
especially  in  the  case  of  cirriforn  clouds.  For  cirrus  clouds 

-3 

with  an  ice  concentration  of  0.01  g*m  a thickness  of  al- 
most 5 km  is  required  to  assure  that  the  assumption  of 
opaqueness  is  reasonably  satisfied.  Thin  trailing  edges  of 
cirrus  which  generally  exhibit  lov;  emissivity  may  be  classi- 
fied, for  example,  in  the  700-400  mb  layer  rather  than  the 
layer  above  400  nb.  For  this  reason,  it  is  suggested  by 
Koffler  et  al.  ['hi]  that  the  cloud  data  be  interpreted  as 
lying  above  700  mb  rather  than  exactly  within  the  700-400  mb 
layer.  If  assumption  number  (4)  is  not  valid  over  a given 
NMC  grid  area  (12  km  at  45°N)  , then  there  will  be  errors  in 
the  calculation  of  the  ar.ount  of  low  cloud.  To  insure  that 
cart))  surface  view  is  not  included  as  low  clouds,  a factor  of 
S^C  is  subtracted  from  the  surface  temperature  T^^  supplied  by 
NMC . 

Results  obtained  by  Koffler  et  aL  [73  ] show  good  large- 
scale  agreem.ent  with  surface  observations  and  pilot  reports 
while  offering  the  advantage  of  greater  geographical  cover- 
age. The  authors  state  that  "problems  and  errors  of  the 


dot 
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system  are  probably  no  greater  than  that  of  the  surface 


observation  system."  There  is  however  a tendency  to  over- 
estim.ate  cloud  amount  and  underestimate  tho*  temperature  of 
cloud  top  lieig);ts  as  a result  of  satellite  sensor  resolution. 
If  the  size  of  a cloud  elcmc.nt  is  near  or  below  the  sensor 
resolutio.n,  the  IR  temperature  spot  will  represent  a com- 
bination of  the  temperature  of  the  surface  and  the  tempera- 
ture of  the  cloud,  thus  causing  the  entire  area  to  be  classi- 
fied as  cloud  v.’ityi  the  cloud  top  height  placed  at  lower  than 
actual  level.  This  type  of  error  can  bo  expected  over  areas 
of  fair  v.'cather  cumulus  typical  of  oceanic  regions.  The 
given  clasr>if ication  system,  can  be  applied  over  oceanic  re- 
gions in  the  Southern  Rcr.isphore  whore  IhhC  data  is  not  avail- 
able by  using  vertical  tcs.peraturc  profile.s  detcrm.incd  from 
Eca-surface  temp- •ratu'-c  csrim.r.tcs  gc.neratcd  from  satellite 
IR  data  as  dorcriyed  in  Rao  [to]. 

historically  or.-,  ef  the  first  attempts  at  m.ultispcctral 
cloud-typo  iCentif icaticn  was  the  work  of  Greaves  a: d Chang 
[ 7C]  v.-’nlcy-i  in"ludcd 


(1) 


an  an  .lys-if:  C'f  tl'.e  statistical  distributions  of 
fiva  cloud  typos  using  thrcc-channel  radiometric 
data  at  sg.’cctral  intervals  .2-4.0  pm,  G.4-6.9  pm, 


and  10-11  pm., 

(2)  derivation  of  cloud  signatures  from  mean  and  stand- 
ard deviations  observed  in  the  reflectance  cliannel 
(.2-4.0  pm.)  and  the  temperature  channel  (10-11  pm), 

(3)  imp  lem.cn  tat  ion  and  testing  of  a decision  system 
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based  on  rectangular  decision  regions,  and 
(4)  formal  development  of  a decision  matrix  system  for 
elliptical  decision  regions. 

Data  from  Nimbus  II  at  pieces  vdicre  a relatively  uniform 
cloud  cover  v:as  reported  over  a large  area,  where  there  v;as 
no  rapid  clearing  or  buildup  of  clouds,  and  where  ground- 
observed  documentation  of  cloud  type  and  cloud  cover  was 
available  vrero  divided  into  a tost  set  and  a desicn  set 


which  v;as  used  for  development  of  cloud -type  signatures  for 


five  different  cloud  types  categorized  as  follov.’s: 

(1)  c'umulus 

(2)  stratus  and/or  stratocumulus 

(3)  altocumulus 

(4)  cirrus 

(5)  cumulonimbus. 

An  analy.'-.is  of  the  statistical  distributions  of  each 
cloud  type  in  each  of  the  throe  radiometric  channels  revealed 
that  while  the  distributions  in  the  rcfiectanco  and  tempera- 
ture c’riannel:;  exliibitod  characteristic  r.axim.a  and  spreads  per 
clcu.d  type,  th.e  information  content  of  distributions  in  the 
water  vapor  channel  (o.4-G.9  pm.)  was  low  in  term.s  of  cloud-type 
differentiation.  A primary  advantage  of  incorporating  tem- 
perature data  from  the  water -vapor  channel  into  a given 
classification  schcT.'.o  was  the  fact  that  cloud  data  could 
easily  be  separated  from  clear  skies  by  classifying  all 
points  as  clear  in  w)iich  the  tcm.ncraturc  obtained  from  the 
watcr-vapor  channel  wa.s  greater  than  240''K.  This  tlircshold 
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proved  adequate  even  in  identification  of  tenuous  cirrus  as 
cloud  data.  VJith  the  assumption  that  the  distributions  are 
approxiraatcly  Gaussian,  moan  and  standard  deviations  for  the 
five  cloud  types  in  the  reflectance  and  temperature  channels 
wore  computed  as  shov.’n  in  Table  2.2  below. 


Cloud  Type 

T 

0 (T) 

R 

o(R) 

cumulus 

274°K 

CO 

o 

0.21 

0.06 

stratus/stra- 

tOCUITiUluS 

281°K 

CO 

o 

0.59 

0.06 

altocumulus 

262'^K 

5.0°K 

0.54 

0.08 

cirrus 

295°K 

3.6°K 

0.10 

0.03 

cum.ulonimbus 

224°K 

6 .4°K 

0.62 
1 

0.09 

Table  2.2:  Cloud-type  signatures  used  by  Greaves  and 

Chang  [ 30]  • 

T represents  mean  temperature  in  the  spectral 
interval  10-11  vm 

0 (T)  represents  standard  deviation  in  the  spec- 
tral interval  10-11  pm 

K represents  m*oan  reflectance  in  the  spectral 
interval  .2-4.0  pm 

0 (R)  represents  standard  deviation  in  the  spec- 
tral interval  .2-4.0  pm 
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It  should  be  noted  that  the  relatively  high  value  of  T for 
cirrus  is  a rc:;sult  of  transmission  to  the  satellite  of  ra- 
diation from  lower  levels.  The  standard  deviation  o(T)  which 
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represents  the  degree  of  nonuniformity  in  cloud-top  height  or 
cloud  thickness  is  seen  to  be  maxiraal  for  cumulonimbus  obser- 


vations . 

A given  data  point  with  reflectance  R and  temperature  T 
was  classified  as  belonging  to  cloud  type  i if 
|R  - R^l  i K 0 (R^) 

and 

It  - T.  I ^ i:  o (T. ) 

where  , T^,  a (R^) , a(T^)  are  values  given  in  Table  2 above 
and  K is  a range  factor  chosen  as  K = 2.5  for  the  given  imple- 
mentation. The  authors  suggest  that  the  optimum  value  of  K 
for  classification  of  most  data  points  v/ithout  excessive  over- 
lap bctv.’oen  cloud  types  lies  between  2.0  and  2.5.  In  cases 
where  a data  point  fell  into  overlapping  regions  of  two  diff- 
erent cloud  typos,  it  was  considered  to  be  correctly  classi- 
fied if  cither  of  the  cloud  types  corresponded  to  the  ground- 
truth  identification  (v.'hich  was  based  on  data  from  ground 
observations  as  well  as  analysis  of  satellite  photographs) . 

The  percentage  of  correctly  classified  data  points  per 
cloud  type  were  as  follows: 

(1)  431  for  cumulus 

(2)  78?.  for  stratus/stratocumulus 

(3)  12%  for  altocumulus 

(4)  B0%  for  cirrus 

(5)  34%  for  cumulonimbus 

The  authors  attribute  the  lack  of  success  in  point-by-point 
identification  of  cumulus  and  cumulonimbus  to  the  nonhomogen- 
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eous  nature  of  those  cloud  fields  and  attribute  the  high  de- 
gree of  success  in  recognition  of  the  remaining  cloud  types 
partially  to  the  fact  that  the  test  samples  consisted  only  of 
regions  identified  as  being  covered  by  a single  type  of  cloud- 
iness. They  suggest  that  spatial  frequencies  be  considered  to 
facilitate  the  identification  of  cuiriuliform  cloud  types. 

A decision  m.atrix  approach  to  cloud  classification  using 
elliptical  boundaries  was  formulated  for  a 2-group  cloud-type 
model  as  follows:  Store  in  the  computer  memory  a 20x20  de- 

cision matrix  D where  the  value  of  the  decision  matrix  at 
D(r,t) [where  r is  the  rth  reflectance  interval  (with  each  of 
20  intervals  representing  a 5%  increment  in  reflectance)  and 
t is  the  t^'^^  temperature  interval  (v;ith  each  of  20  intervals 
representing  a 5%  increm.ent  in  temperature)]  can  assume  any  of 
the  values  0,  1,  2,  3 which  stand  for  the  four  possible  classi- 
f ications 

0 - no  decision 

1 - cloud  type  1 

2 - cloud  type  2 

3 - ambiguous  between  cloud  types  1 and  2 

The  value  assigned  to  clement  D(r,t)  of  matrix  D depends  on 
whether  or  not 


R-R.  ~ 

2 

I — 

H 
I 

X 

+ 

— 

•H 

D 

I 

o (T^) 

2 

where  N is  a range  factor  (sim.ilar  to  the  factor  K above)  , 
R = 5r 
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T = 5t  + 200° 


^ represent  cloud-type  signatures 

similar  to  those  in  Table  2 above. 

If  the  above  inequality  is  satisfied  for  both  group  1 (i=l)  and 
group  2 (i=2)  , then  the  value  of  D(r,t)  is  3 v;hich  means  that 
the  classification  is  aimbiguous.  Once  the  decision  matrix  is 
stored  in  the  computer,  the  classification  of  any  data  point 
is  achieved  by  computing  r and  t from  the  reflectance  R and 
temperature  T and  then  consulting  the  entry  D(r,t)  in  the  de- 
cision matrix  for  the  corresponding  classification.  The  authors 
suggest  the  decision  matrix  approach  for  future  configurations 
for  cloud-type  classification. 

Shenk  and  Holub  [12]  proposed  a multispectral  cloud-type 
identification  method  based  on  a cloud-type  decision  matrix 
j for  the  four  spectral  regions  of  0. 2-4.0  ym,  6.4-6. 9 ym,  10- 

11  ym,  and  20-23  ym.  Decision  categories  v.’ere 


1 
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(1)  cirrus  clouds 

(2)  cirrus  with  lower  clouds 

(3)  cumulonimbus  and/or  cirrostratus 

(4)  cum.ulonimbus 

(5)  middle  clouds 

(6)  middle  clouds  with  cirrus  above 

(7)  stratus  or  stratocumulus 

(8)  cumulus 

(9)  clear 

(10)  no  decision 

A data  point  from  Nimbus  3 satellite  imagery  was  classified  in- 
to one  of  the  categories  1-9  if  measurements  in  all  four  spec- 
tral channel  fell  between  the  lower  and  upper  threshold  for  that 
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particular  cloud  type.  Otherwise,  the  data  point  fell  into  the 
no  decision  category.  The  authors  recommended  examining  adjacent 
cloud  decisions  to  obtain  information  on  cloud  conditions  at 
"no  decision"  points. 

Threshold  values  were  calculated  from  a radiative  trans- 
fer model  which  accepted  as  inputs  average  temperature  and 
moisture  sounds  for  the  given  test  area  during  the  study  period, 
a climatological  ozone  profile,  and  cloud  tops  at  different 
heights.  It  was  assumed  that  lov,-cloud  tops  were  at  750  mb, 
middle  cloud  tops  at  450  mb,  cirrus -cloud  tops  above  300  mb, 
and  cumulonimbius -cloud  tops  below  225  mb. 

From  a total  of  611  cloud -type  decisions  over  a 30-mile 
swath  covered  by  aircraft  photographs,  85  decisions  were  mis- 
classified.  Two-thirds  of  the  misclassif ications  represented 

a confusion  between  groups  1 and  2,  which  probably  could  have 
been  avoided  by  different  threshold  choices  for  the  visual 
channel.  The  remainder  of  the  misclassif ications  except  for 
two  points  were  the  result  of  labelling  scattered  low  clouds 
or  clear  conditions  as  cirrus.  Some  of  the  misclassif ication 
errors  were  attributed  to  earth  location  errors  in  the 
satellite  data,  spatial  resolutions,  and  time  difference  be- 
tween aircraft  and  satellite  observations.  Modifications  in 
cloud-type  categories  and  cloud -type  thresholds  were  later 
adopted  (Shenk  et  al.  [ 7Q)  which  significantly  in\proved 
classification  accuracy. 

Coincident  infrared  and  visual  data  from  NOAA-1  satellite 
were  used  by  Booth  [73]  to  test  the  discriminatory  ability  of 
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Bayesian  and  maximiun-likelihood  classifiers  to  separate  2100 
cloud  observations  into  a 6-group  tropical  cloud -type  model 
consisting  of 

(1)  clear  skies 

(2)  cumulus 

(3)  stratocumulus 

(4)  cumulonimbus 

<5)  cirrus 

(6)  cirrus  with  lower  clouds 

and  into  a 5-group  tropical  cloud-type  model  obtained  from  the 
6-group  model  by  com±>ining  cum.ulus  and  stratocumulus  into  one 
group.  A comparative  study  of  the  effect  on  each  of  the 
classifiers  of  selecting  feature  vectors  from 

(1)  single-channel  visual  (.5-. 7 ym)  data 

(2)  single-channel  infrared  (10.5-12.5  ym)  data 

(3)  dual-channel  visual  and  infrared  data 
revealed  that  accuracy  using  single-channel  infrared  data 
approached  that  obtained  by  dual-channel  classification. 

A combination  of  features  representing  32  spatial  distri- 
bution measurements  and  32  spatial  frequency  measurements  w'ere 
calculated  for  each  32x32  observation  matrix  (approximately 
54x96  nmi.)  of  coincident  visual  and  infrared  data.  The 
spatial  distribution  measurements  selected  for  both  IR  and 
visual  data  were  mean,  standard  deviation,  coefficient  of 
variation,  skewness,  kurtosis,  range,  mean-median,  primary 
mode  / secondary  mode,  average  gradient,  maximum  quadrant  co- 
efficient of  variation-minimum  quadrant  coefficient  of  varia- 
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tion,  and  cumulative  frequency  values  at  1%,  16%,  50%,  84%, 


and  99%  levels.  The  power  spectrum  measurements  selected  for 
bot)i  IR  and  visual  data  were  values  at  V7ave  numbers  1-11  and 
linear  combinr.tions  of  various  v,*ave  numl:iers  from  1-11.  The 
spatial  distribution  measurements  have  typically  been  used  to 
describe  size,  shape,  and  texture  of  cloud  images.  The  power 
spectrum  measurements  denote  the  amount  of  variance  in  an 
observation  due  to  clouds  of  various  sizes.  For  example,  given 
a grid  size  of  54x96  nmi . , wave  number  1 corresponds  to  a cloud 
size  of  27x48  nmi. 

Both  the  Bayesian  and  maxim.um-likelihood  classifiers  v;ere 
designed  on  feature  vectors  obtained  from  1050  observation  ma- 
trices (half  of  the  sample  set)  and  were  tested  on  both  the 
design  set  (dependent  data)  and  the  rem.aining  half  of  the  sam- 
ple set  (independent  data)  not  used  for  design.  The  con- 
ditional probabilities  p(i|x)  where  p(il>:)  is  the  posteriori 
probability  that  an  observation  belongs  to  group  i given 
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observation  vector  X wore  calculated  under  the  assumptions 


that  the  feature  vectors  wore  normally  distributed  and  that 
the  covariance  matrix  of  each  class  could  be  approximated  by 


the  pooled  within-groups 


dispersion  matrix  ) where 


where  ] . 


is 


(Fi-DljH-  (F^-l)^^-f...+  (F^.-l)I^^ 
N-k 

the  covariance  matrix  for  group 


i , 


Fj^  is  the  group 


relative  frequencies  obtained  from  the  entire  sample  set  (in- 
dependent and  dependent  data) . For  the  maximum-likelihood 
classifier  all  calculations  were  performed  as  for  the  Bayesian 
clas:;ifier  with  the  exception  that  equal  a priori  probabilities 
were  used.  The  decision  rule  was  to  classify  and  observation 
into  group  i if 


P(i|X)  > P(jl>:)  for  all  j 7^  i 


where 


P(i|x) 

I P(xli)  -P(i) 
i=l 


P(xli)  exp[-  |(X-I.:  )'^’L ^(X-M.)] 


IJpl  dofi.ned  as  t}jc  determinant  of  the  pooled  dispersion 
matrix 

k defined  as  the  number  of  classification  grovips 

M.  defined  as  the  iv,can  vector  for  the  ith  classifica- 
1 

t,ion  group 

n defined  as  the  niunber  of  features  in  feature  vector 


With  the  threshold  level  for  screening  set  at  the  1% 
significance  level,  thirty-two  features-  (sixteen  from  each 
channel)  for  dual-channel  data,  sixteen  for  single-channel 
IR  data,  and  seventeen  for  single-channel  visual  data  entered 
into  the  discriminant  functions  g(X)^  (where  P(x|i)  = exp[g(X)p]) 
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calculated  for  the  maximum-likelihood  and  Bayes  classifiers 
using  the  stepv.’ise  multiple-discriminant  analysis  program 
(EMDOVM)  from  the  UCLA  biomedical  package.  The  feature  v;hich 
was  t};e  best  single  discriminator  for  dual-channel  data  v;as 
the  IR  viilue  at  tl'.e  11  cumulative -frequency  level  (at  the 
cold  cr.d  of  the  cumulative- frequency  curve)  v;hich  represents 
(for  observations  containing  clouds)  the  temperature  of  cloud 
tops.  For  opaque  clouds  the  thirty-two  features  chosen  tend 
to  measure  various  characteristics  of  cloud  tops;  however,  for 
semi-transparent  clouds  such  as  cirrifonr  the  effects  of 
underlying  visual  and  infrared  radiation  contaminate  the 
measurement  of  cloud- top  characteristics.  The  miinimum  mean 
temperature  ranges  were  found  in  stratccumulus  samples.  The 
m.aximium  mean -temperature  ranges  were  found  in  cumiulonimibus  and 
cirrus  observations  with  the  coldest  mean  temperature  found  in 
cum.ulonim.bus  observations.  T)ie  moan  tem.porature  for  cirrus 
observations  was  relatively  high  as  a consequence  of  radiation 
transmitted  from  below  the  cloud  tops.  The  m.aximumi  mean  tem- 
perature gradient  v.-as  found  in  cirrus  observations.  The 
feature  whicJi  v;as  the  second  best  discriminator  for  dual- 
channel  data  was  ti',e  visual  value  at  the  991  cumulative  fre- 
quency level  vdiicli  represents  high  albedo.  Kean  relative 
brigl'.tness  was  a naxi.mum  for  cumulo.nimljus  observations , a 

or  cumulus 
ature  of 
tratocum.ulus 
s maximum  for 


cumulonimbus  and  minimum  for  observations  containing  only 
cirriform  clouds.  The  largest  mean  briglitness  gradients  were 
found  in  cumuloninibus  and  stratocumulus  observations. 
Measurements  of  pov/er-sj^cefrum  features  revealed  that  the 
first  five  wave  numljers  accounted  for  791  of  the  total 
variance  in  cumuliform  observations  and  89i  of  the  total 
variance  in  cumulonimbus  observations. 

Difference  in  overall  performance  bctv.’cen  the  maximum- 
likelihood  classifier  and  the  Bayes  classifier  did  not  vary 
by  more  than  3%  no  matter  whether  dependent  or  independent 
data  v;as  being  considered  or  whether  single-channel  IR,  single- 
channel  visual,  or  dual-channel  data  was  being  considered  or 
v.’hethcr  a 5-croup  or  6-group  cloud- type  model  was  considered. 
For  dual-channel  Bayes  classification  for  the  6-group  model 
the  overall  hit  percentage  was  76-1  for  the  dependent  (design) 
data  and  631  for  the  independent  data;  for  single-channel  IR 
Bayes  classified’. tion  for  the  6-group  model  the  overall  hit 
percentage  was  601  for  dependent  data  and  571  for  independent 
data;  for  single-channel  visual  Bayes  classification  for  the 
6-group  model  tlic  overall  hit  percentage  was  541  for  dependent 
data  and  431  for  independent  data.  In  general,  classification 
results  v.'erc  better  using  the  Bayes  classifier  on  groups  with 
high  a priori  probabilities  whereas  the  overall  hit  percentage 
lor  groups  with  lov.’  a priori  probabilities  was  better  for  the 
ma.ximum-li):eli):ood  classifier  than  for  the  Bayes  classifier. 

The  major  misclassif ications  occurred  in  distinguishing 
cunulrnir.L'Us  from  cirrus  with  lower  clouds,  cirrus  from  cumu- 
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lonimbus,  cirrus  with  lov;or  clouds  from  cirrus,  and  camulus 
from  clear  skies.  Relatively  few  curaulonimbus , cirrus,  and 
cirrus  with  lower  clouds  were  misclassif ied  as  clear  or  into 
one  of  the  low-cloud  groups.  The  pattern  of  errors  obtained 
might  be  expected  as  a result  of  the  procedure  for  manual 
classification  of  cloud  typos  by  the  trained  meteorologists. 
If  an  observation  contained  cumulonimljUG  clouds  and  cirriform 
clouds,  it  was  labelled  as  cumulonimijus  regardless  of  cloud 
amount  of  each  of  the  two  types.  Also  vdien  cumulonimbus  and 
cumulus  clouds  occurred  in  the  same  observation,  it  was 
labelled  as  cumulonimbus,  and  when  an  observation  contained  a 
small  eimount  of  cumulus  cloud,  it  was  labelled  as  cumulus 
rather  than  clear  skies. 

The  overall  hit  percentages  for  the  Bayes  classifier 
applied  to  features  extracted  from,  single-channel,  IR  data, 
single-channel  visual  data,  and  dual-channel  data  are  given 
in  Table  2.3  below. 


Six -Group  Model 

Five-Group  Model 

Channel 

Dependent  Data 

Independent  Data 

Dependent  Data 

Independent  j 

Data  ! 

Visual 

5^.?, 

43% 

55% 

44% 

IR 

60% 

57% 

69% 

63% 

Visual/IR 

76% 

63% 

81% 

69% 

Table  2.3:  Overall  Hit  Percentages  of  Bayes  Classifier  obtained  by 
Booth  [73]. 
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Thus  it  can  be  seen  that  features  from  the  single-channel  IR 
data  led  to  overall  better  classification  results  than  could 


be  obtained  using  features  only  from  single-channel  visual 
data.  Hov;ever  clear  skies  can  best  bo  identified  using 
single-channel  visual  data.  Cumulus  clcuds  V/cre  best  identi- 
fied by  single-channel  IR  data.  For  accurate  identification 
of  cirrus  clouds  dual-channel  data  performed  best. 

Recommendations  made  by  Booth  for  future  research  in 
cloud  classification  techniques  included: 

(1)  isolation  of  cloud  elements  in  an  observation  by  use 
of  sea-surface  temperature  measurements  vith  sub- 
sequent calculation  of  features  only  on  regions  of 
an  observation  vhich  contain  cloud  elements  rather 
than  on  the  entire  observation 

(2)  addition  of  a "no-decision"  class  j.n  which  those 
observations  whose  a posteriori  probability  was  not 
greater  than  a given  threshold  could  be  placed 

(3)  additional  training  sam.ples  for  the  groups  of  strato- 
cumuluG  and  cirrus  with  lov.’cr  clouds 

(4)  investigation  of  non-parametric  methods  of  discrimi- 
nant analysis. 
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3.2.  Cl oud  Sogmen tci_t u'ji 

Scene  scymontation  refers  either  to  the  partial  or  complete  decomposi- 
tion of  image  data  into  specific  parts.  The  simplest  example  of  segmentation 
of  meteorological  satellite  data  is  the  separation  of  cloud  from  no  cloud 
points,  hartitioning  image  data  into  cloud-type  regions  is  another  example 
of  segmentation.  Partial  segmentation  of  an  image  may  result  in  the  extrac- 
tion of  cloud  objects  or  patterns  such  as  thunderstorm  cells,  low-level 
cumulus  cloud  line  feeder  bands,  or  arc-shaped  mesohigh  boundaries. 

Basic  segmentation  techniques  for  decomposition  or  extraction  of  objects 
include  thresholding  local  properties,  edge  detection,  template  matching, 
raster  tracking,  and  region  growing.  Segmentation  of  meteorological  satel- 
lite data,  either  for  cloud- type  or  cloud-pattern  extraction  or  for  observa- 
tion of  thunderstorm  structure,  has,  for  the  most  part,  been  the  result  of 
application  of  thresholding  or  edge  detection  techniques.  Matching  tech- 
niques have  been  used  to  track  the  movement  of  cloud  patterns.  Raster  track- 
ing and  its  variations  (such  as  omnidirectional  tracking,  etc.)  are  employed 
in  algorithms  for  line  and/or  curve  detection.  Region  growing  techniques 
start  with  initial  partitions  of  an  image  into  homogeneous  regions  (for  ex- 
ample, regions  of  constant  gray  level)  and  then  meme  or  split  regions  on  the 
basis  of  so.ne  criterion  for  the  "goodness"  of  a partition.  Merging  of  par- 
titions may  continue  until  the  gray  level  variance  exceeds  a specific  thres- 
hold or  may  be  based  on  more  complicated  semantic  criteria  (see  Yakimowsky 
and  Feldman  [1973]). 

The  selection  of  appropriate  cloud  pattern  thresholds  is  one  of  the  most 
important  areas  of  investigation  for  extraction  of  severe  storm  features. 
Parameters  defining  the  vertical  growth  rate  of  a cell,  sneed-of -cel  1 move- 
ment, presence  of  merging  cells,  and  presence  of  splitting  cells  depend  criti- 
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cally  on  appfopi'ioLo  tC'inperature  coiUourinq  of  limited-scan  infrared  data. 
According  to  Adler  [^576],  "t‘ie  area  expansion  inside  particular  isotherms  is 
a more  accurate  and  more  sensitive  parameter  [of  thundcrstonii  growth]  in  all 
but  the  very  early  stages"  Uian  tlic  minimum  hlackliody  temperature.  Adler 
deaionstrated  that  small  areas  of  colder  temperatures  within  the  anvil  often 
exhibited  fluctuations  associated  with  convective  activity  that  were  not  ob- 
served in  the  expansion  of  tiic  cirrus  anvil. 

Automatic  selection  of  (iray-lcvel  thresholds  is  usually  approached  by 
assuming  that  cbjecLs  or  patterns  in  an  image  can  be  characteri ?cd  by  modes 
of  normal  distributions  in  gray-level  histograms.  For  the  mode  method,  the 
threshold  is  selectc'c:  at  the  bottom  of  a valley  between  two  peaks.  If  the 
valley  is  broad  and  the  pc.;ks  are  of  unequal  size,  a Laplacian  ooerator  can 
be  apoiijci  to  the  nirture  to  determi r'-*  y-.e  gt'ay-level  value  of  points  that  lie 
on  or  n‘’Tr  object  Lm>rders.  Automatic  threshold  selection  of  cloud  pictures 
(visiPle  cl'’ta;  ur.m:n  a Margua.rdt  rearcssicn  procedure  was  investigated  by 
Yen  [K'7-'].  The  procedure  assun;ed  tliat  the  gray-level  distribution  of  the  pic- 
ture was  a sum  oi  normals.  Tlic  success  of  the  method  was  strongly  dependent 
on  I'u'  number  of  i;ormal  density  functions  (two  or  three)  postulated.  Object/ 
background  threshold'',  were  selected  at  the  intersections  of  the  first  and  se- 
cor.d  nonaal  densiiv  functions. 

kdge  tint oc*  ioc.  oooratc'rs  have  also  been  applied  to  mcteorologi cal  satel- 
lite d.ila  to  extract  clcucls  from  background.  Smootliinq  or  averaging  the  image 
prior  to  appl  ica  ticwi  of  edgu'  detection  operators  can  soiiictimes  improve  the 
results.  The  Fobci'ts  digital  gradient  operator,  which  is  defined  for  any 
2x2  array  of  gr.''y  level  values  ^ ^ as  max  ( ! A-P | , j B-C  ] ) , often  performs 
as  well  as  more  e>;i''‘irsi ve  or  sophisticated  edge  techniques  on  specific  cloud 
images.  Edges  can  also  be  detecter'  by  high-pass  spatial  frequency  filtering. 
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The  output  of  cdgo  detection  operators  can  bo  used  in  conjunction  with  border 
tracking  algorithms  to  outline  an  object  or  the  output  can  be  used  to  deter- 
mine a gray-level  tlireshold  (isotiierm)  for  contouring  an  object. 

Segmentation  of  infrared  data  into  cloud-type  objects  by  analysis  of  the 
value  of  ti'O  average  border  edge  strengtii  of  connected  components  at  each 
possible  tCM'.perature  value  is  discussed  in  i’arikh  [l977b] . Three  edge  operators 
were  developed.  tdc;e  o'serators  v/hich  differentiated  between  stop  edges  and 
ramp  edges  were  more  sensitive  to  isolation  of  uni  form- textured  cloud  patterns 
than  operators  which  searched  only  for  maximum  difference  between  neighbors. 
Results  of  autoisatic  throsiiold  selection  based  on  this  hybridization  of  thres- 
hold and  edeo-detection  techniques  could  probably  be  improved  by  procedures 
such  as  using  a weighted  average  border-edge  strcngtfi  feature  with  less  weight 
at  border  points  wi  tli  fuzzy  edges. 

Techniques  for  automatic  selection  of  temperature  contours  for  analysis 
of  thunders toiT’  structure  need  to  be  developed.  These  techniques  could  be 
used  to  prccTico  enhanced  infrared  imagery  and  to  extract  features  which  re- 
late to  severe  thunderstorm  activity.  The  selection  of  appropriate  tempera- 
ture contours  could  bo  based  on  analysis  of  features  otlier  than  border-edge 
strength.  Features  such  as  standard  deviation  of  temperature  (within  tlie  con- 
toured areas),  uniformity  of  texture,  moan-visible  brichtness,  etc.,  should 
also  be  examined. 

Streaks,  curves,  and  lines  can  be  detected  by  examining  local  neiglabor- 
hoods  (points  or  adjacent  nonovcrlappi ng  averages  of  points)  for  line-like  or 
streak-like  patterns.  Results  of  streak  detection  operations  can  often  be 
significantly  improved  by  repeating  tfic  operation  more  than  once.  A descrip- 
tion of  linear  and  nonlinear  operators  for  streak  and  line  detection  can  be 
found  in  Rosenfold  and  Kak  0976].  Similar  operators  can  be  designed  to  de- 
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tect  intersections  or  crossings  of  cloud  patterns.  A comparative  study  of 
curve  detection  operators,  line  detection  operators,  and  operators  for  detec- 
tion of  branching  or  crossings  of  curves  needs  to  be  conducted  in  order  to 
determine  the  most  effective  pattern  recognition  techniques  for  identifica- 
tion of  mesohigh  boundary  interactions  and  low-level  cumulus  cloud-line  feeder 
bands . 


3.3.  C]  0 u d_G c 0 me  try 

Measures  for  describing  geometrical  properties  or  relationships  of  cloud 
cells  form  an  integral  part  of  pattern  recognition  algorithms  for  detection 
of  severe  storm  features  such  as  merging  and  for  splitting  cells,  comma-shaped 
patterns,  oval-shaped  leading  edges  of  mesohighs,  and  vertical  growth  rate 
of  cells.  Distance  measures,  measures  of  elongation,  curvature  measures  and 
area  measures  must  be  specified. 

The  concept  of  merging  and/or  splitting  cells  implies  that  a distance 
measure  between  two  cells  exist.  The  distance  between  two  cells  can  be  a 
function,  for  example,  of  the  distance  between  the  centroids  of  the  cells, 
the  minimum  distance  between  edge  points,  or  the  distance  between  the  mini- 
mum temperatures  of  the  cells.  If  merging  of  two  cells  is  defined  in  terms 
of  a zero  distance  between  one  or  more  edge  points  of  the  cells,  then  the 
recognition  of  merging  and/or  splitting  cells  would  be  dependent  on  the  selec- 
tion of  an  appropriate  isotherm  to  contour  the  cells. 

Techniques  for  location  of  comma-shaped  patterns  (such  as  thinning,  thick 
curve  following,  etc.)  could  bo  developed  or,  alternatively,  measures  could 
be  defined  to  detect  comma-like  shapes.  A possible  definition  of  a comma 
pattern  might  involve  location  of  ttic  medial  axis  of  the  pattern  and  examina- 
tion of  ratios  of  the  widths  to  the  length  of  the  pattern. 


r 

Ovol-shapcd  boundaries  of  niesohiqhs  could  be  defined  in  terms  of  arcs 
v/ith  expected  curvature  values  greater  than  a given  thresliold.  The  curva- 
ture of  a point  on  an  arc  is  given  by  the  difference  between  the  left  and 
right  slopes  at  the  point.  For  an  arc  consisting  of  points  po,Pi  , . . . ,Pp > 
the  left  and  right  k-slopes  at  the  point  Pi  are  the  directions  from  p.j  to 
Pi_k  and  from  p^  to  p-j+|<,  respectively. 

I The  area  of  a digital  object  is  defined  as  the  number  of  points  in  the 

I 

' object.  In  order  to  calculate  che  vertical  growth  rate  of  cells  in  terms  of 

[ 

I a cloud  expansion  parameter,  appropriate  isotherms  for  cell  definition  must 

be  selected  throughout  the  time  period  of  observation.  The  selection  of 
the  warmest  temperature  that  is  unambiguously  part  of  the  entity  is  discussed 
in  Adler  [1976J. 

3.4.  Cloud  Tracking 

The  extrapolation  of  time-dependent  features  is  a fundamental  require- 
ment for  prediction  of  severe  weather.  Algorithms  for  automatic  selection 
I and  tracking  of  cloud  objects  have  been  developed  primarily  for  wind  velocity 

I estimation.  These  algorithms  are  reviewed  in  Sections  3.4.1  and  3.4.2.  Auto- 

I mated  algorithms  for  tracking  thunderstorm  cells  and  other  similar  time- 

I varying  features  remain  to  be  developed.  The  time-varying  nature  of  these  fea- 

tures distinguishes  the  problem  of  recognition  of  severe  storm  features  from 
most  current  problems  in  pattern  recognition  and  imaae  processing. 


..A,. 


3 . 4 . 1 . Automatic  Selection  Techniques  for  Cloud  Targets 

Automatic  methods  for  wind  velocity  estimation  which  do 
not  rely  upon  the  skill  of  a trained  meteorologist  to  select 
cloud  targets  to  ))c  tr<;ickc?d  from  frame  to  frame  tend  to  imple- 
ment a particular  aspect  of  tlic  manual  cloud-target  selection 
procedures  Prior  to  mamial  selection  of  a particular  indivi- 
dual elcniicnt  or  feature  to  bo  traced,  a meteorologist  will 
usually  a'jsess  the  synoptic  situation  and  mentally  note  the 
general  cloud  pattern  flow  over  a largo  area.  Automatic 
cross-correlation  metiiods  integrate  cloud-pattern  motion  over 
an  area  (usually  32x32  or  64x0^  picture  points).  The  binary 
matching  technique  developed  by  Bristor  [72]  emulated  the 
manual  procedure  of  tracking  cloud  edges  from  fram.e  to  frame. 

As  a result  of  factors  sucli  as  evaporation,  not  all  edges  are 
conserved  from  frame  to  frame.  In  practice  many  meteorologists 
would  choose  to  follow  brigjitness  centers  (not  geometric  cen- 
ters) of  cloud  masses  from  frame  to  frame  rather  than  follow 
edges.  Tlie  clustering  algorithm  of  Endlich  et  al.  [71]  in 
w’liicli  cloud  data  were  represented  by  brightness  centers  is  an 
analog  to  this  procedure.  Lo  et  al.  [74]  emphasized  Fourier 
transform  filtering  techniques  for  selecting  individual  cloud 
elements  of  the  same  approximate  size  (preferably  about  10  Iuti) 
as  those  normally  used  in  manual  selection  tecliniques. 

The  operational  model  of  the  National  Environmental  Satel- 
lite Service  (NESS)  described  by  liradford  et  al,  [72]  » computes 
cloud-displacement  vectors  from  64x64  cross-correlation  ma- 
trices resulting  from  Irackward  cross-correlation  of  32x32 
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cloud- target  areas  from  the  second  frame  with  G4x64  areas  from 
the  first  frame.  Ecich  32x32  array  (and  each  64xG4  array)  is 
centered  at  a five-degree  latitude-longitude  intersection. 

Tile  ciioict;  of  specific  sizes  for  the  two  arrays  was  based  on 
tile  expected  siieed  oJ'  the  cloiuis.  The  use  of  a smaller  array 
for  the  target  are^l  than  for  the  search  area  was  found  to  im- 
prove the-  accuracy  of  cross-correlation  wind  velocity  estimates 
for  tliosc  areas  in  wiiich  clouds  were  m.oving  in  or  out  of  the 
boundaries  or  in  wh.icii  the  l^oundary  between  two  different  cloud 
patterns  came  inside  the  G4x64  grid. 

Bristor  [’^2]  developed  a i^inary  matching  technique  for 
determining  cloud -d i splacement  vectors  as  early  as  19G7.  A 
cloud  edge  ■.-.'as  defined  by  assigning  to  it  a brightness  isopleth. 
Tlio  brigiitnoss  range  of  the  edge  lay  betv.'een  the  brighter 
values  witliin  a cloud  and  the  darker  values  of  tlie  background 
area.  The  isopleth  may  reprc.sont  actual  cloud  edges  or  it  may 
represent  a boundary  v;hore  the  cloud  field  becomes  broken  into 
elements  which  cannot  be  resolved  by  the  satellite  sensor  sys- 


tem. Each,  brightness  image  v;as  reduced  to  a binary  image  with 
c.loud  ociges  depicted  by  I's  and  darker  or  brighter  areas  by 
O's.  Illustrations  of  typical  binary  images  can  be  found  in 


Levs-'  t t al.  (71  I . The  procedure  for  computing  cloud  displace- 
! , f rom  these  binary  images  vv^ill  be  presented  in  the  next 


' . ct  al,  [71]  reduced  cloud  patterns  in  successive 
images  to  centers  of  briglitncss  wliose  displace- 
f-'  frame  agreed  qualitatively  with  cloud 
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motion  vectors  calculated  manually  from  time-lapse  data  of  the 
same  scene.  An  ISODATA  (Iterative  Self-Organizing  Data 
Analysis)  algorithm,  v;as  used  to  select  the  centers  of  bright- 
ness. ISODATA  is  a clustering  procedure  for  obtainiiig  maxin\um 
li  hell  iiood  osLinatcs.  for  mean  values  (centers)  vjith  splitting 
and  merging  or)tions  for  splitting  clusters  with  too  much  within 
cluster  variability  and  merging  clusters  with  too  little 
betwecn-cluster  variability . 

The  program  functions  as  follows.  Scaling  of  tlio  three 
variables  — x ]:>osition,  y position,  and  brightness  value 
B — is  adjusted  in  accordance  v;ith  w’eighting  factors  appro- 
priate  for  tlie  tlirce-dimensional  distance  {Ax'  + Ay  + AB  } ' 
bct’wecn  input  points.  Scales  chosen  for  x and  y values  ranged 
from.  1 to  120  in  increments  of  4 since  the  original  120:-;120 
a^ray  of  satellite  image  points  v.’as  dofocusea  or  averaged  to 
produce  a 3 0x30  array.  The  scales  chosen  for  the  brightnes.s 
value  ranged  from  0 to  150  in  steps  of  10  with  all  points  v;ith 
brightness  values  less  than  60  discarded  as  background  elements 
Then  the  mean  value  (x,  y,  B)  and  the  root  mean  square  distance 

1 

r = K[o^(x)  + p^(y)  + c^(B)]^ 

where 

p 9 2 

o (x)  , o"  (y)  , o (B)  = variances  of  x,  y,  B,  respectively, 

K = sphere  factor,  cfioscn  as  0.7  for  thi.s  study  to  yield 
10-15  centers  per  grid  area 

were  computed  from  the  values  of  all  input  points.  The  first 
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brightness  center  was  specified  as  the  mean  value  (x,  y,  B) . 
Data  points  were  then  compared  in  a sequential  fashion  to  each 
brightness  center  in  order  to  determine  if  the  distance;  from 
the  d.it.i  point  to  the  J:ir i gl; tness  center  is  Ic.'ss  tJ^an  the  radius 
r (defined  above).  If  tlie  distance  is  less  t)ian  or  equal  to 
tJic  radius  r,  the  data  point  is  assigned  to  that  cluster. 
Otlierwisc,  if  the  distance  from  every  brightness  center  is 
greater  them  the  radius  r,  the  f'oint  is  cliosc;n  as  a new  center. 
After  all  [joints  have  'nc-en  processed,  a now  set  of  brightness 
centers  (iac:an  values  of  clusters)  is  computc'd  for  each  group. 

7it  t!iis  st.ige,  spl  itting  and  lum.ping  options  can  be 
applied.  For  sn.litting,  the  standard  deviations  within  each 
grou;i  for  each  of  t;ie  variables  x,  y,  B are  com.puted  and  the 
groLij)  having  tiie  J.argest  value  for  any  one  dimension  is  ten- 
tati'.’cly  spil  i t through  its  uean  value.  Ne-w  centers  are  compu- 
ted, and  if  the  distance  between  the  two  new  centers  is  greater 
than  some  parameter  C^,  the  tv;o  new  clusters  replace  the 
originai  one.  For  lumpin.g,  if  tw'o  cluster  centers  arc  closer 
th.an  the  two  clusters  are  merged  into  one.  It  was  found 

experi  mental]  y llial  th.e  splitting  anti  merging  operations  did 
not  i.iiprove  the  accuracy  of  cloud  motion  vectors.  The  next 
stage  in  the  I FC'DATA  [jrogram  was  to  discard  all  centers  that 
had  less  than  0,,  members,  where  was  set  equal  to  three  for 
this  application.  The  MOTION  program  matched  brightness  cen- 
ters on  tv.’o  picture's  and  com]juted  distances  betwet'n  matched 
pairs  . 

Lo  ct  al.  I 74]  applied  low-pass  filters,  band-pass  filters, 
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nnd  high-pass  filters  of  various  sizes  to  Fourier  transforms 


of  successive  satellites  pictures  in  order  to  separate  cloud 
conpoiu'.'its  of  different  sizess  and  di'seard  cloud  elements  vdiicli 
are  too  Ivirgc;  to  move  passively  v;ith  the'  v;ind.  In  a limited 
regi.o;.  clouds  o'"  different  sizes  tend  to  occur  at  difforemt 
ii  1 L i tud'-'s  and  cons.eejuenti.y  represent  diffc'rcnt:  motion  vectors. 
Results  from  studies  on  simulated  and  rc-al- cloud  data  revealed 
tiiat  the  filtering  technique  did  not  signif icantly  improve 
v;ir,d  velc^city  vcct.ors  for  multi-layered,  mixed-motion  cloud 
regions.  'i'ht'  fiJteiing  teclu-iieiuc  failed  to  catc'gorize  cloud 
size  by  wavelength.  High-frequency  components  contained  in- 
forr’.<'it-ion  not  only  abo'at  small  clouds  but  also  about  sharp 
edees  of  large  clouds  and  similarly  lov;- f roqucncy  components 
contained  informn.tion  not  only  about  large  cloud  masses  but 
also  <ibout  small  clouds  witli  lov;  contrast  and/or  blurred  edges. 
P i f f (":c  n t.  filters  or  different  approaches  to  the  problem  wore 
su<ig<''st '■(!  for  future  research.  In  particular,  thresholding 
p!  nc> -r.si  to  .sepa.rate  different  cloud  typos,  enlianccmont  pro- 
ct'durcr;  much  ns  two-tone  cl  oud/background  images  to  yield 
rduirpci  cro:;s-cer  re  I a i ion  pe,i!;r.,  and  multispcctral  cloud- 
cl  ass  i f icatior.  inc;t!iods  were  expected  to  improve  wind  velocity 
estimates  for  niulti-laycrcd  fields. 
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3.4.2. 


Automatic  Techniques  for  Trackin')  Cloud  Targets 
Automatic  techniques  for  wind-velocity  estimation  such  as 


cross  correlation,  Euclidean  norm,  and  phase  difference  calcu- 
late a cloud  displacement  vector  from  the  bri.cjhtncss  mcasurc- 
mentt;  of  a given  r.ector.  Foi"  these  tcchnicjues , the  entire  sec- 
tor is  chosen  as  the  cloud  target  to  be  tracked  from  frame  to 
frame.  Cross-correlation  methods  can  be  modified  by  selection 
of  a subsc't  of  the  given  sector  to  be  correlated  with  a sector 
from  the  second  image  of  tlie  same  size  as  the  given  sector. 
Ptiasc  difference  and  cross-correlation  Fourier -methods  can  be 
modified  by  filtering  a Fourier  transform  prior  to  calculation 
of  a cloud-displacement  vector.  Binary  matching  techniques 
and  itc:rative  fitting  tecliniques  compute  a cloud  displacement 
vector  from,  cloud  features  or  cloud  targets  derived  from  the 
briglitness  m.casurcmen ts  of  the  given  sector,  i.c.,  from  edge 
features  and  brightness  centers,  respectively. 

For  cross- correlation  techniques  a cloud  displacement 
vector,  represented  by  a displacement  index  of  p picture 
elcm.ents  in  the  x-dircction  (direction  5f  ascending  columns) 
and  a di splaccment  index  of  q picture  elements  in  the  y direc- 
tion (direction  of  ascending  rov;s)  , is  defined  by  the  pair  of 
indices  (p',  q')  fox*  which  tlie  correlation  coefficient 
R(p',  q')  = max[R(p,q)]  for  all  p,  q where  -P  ^ p ^ P and 
-Q  i-  q X.  Q.  Tlie  values  of  P and  Q are  determined  from  con- 
sideiration  of  reasonable  values  for  cloud  displacements . If 
F and  G are  the  means  respectively  of  a sector  f(i,j),  for 
i=l,...,n^  and  j=l,...,nj  from  the  first  image  and  a sector 
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for  and  j = l,...,nj  from  tho  second  image, 

then  tho  correlation  coefficient  R(p,q)  is  defined  as 


v;licre 


R (p,q) 


Cov (p ,q) 

0,0 

f g 


Cov (p ,q) 


1 

{ (n^-p) (n j-q) } -1 


n , -q  n . -p  _ _ 

J^If  (i+p,  j+q)-F]  Ig  {i,  j)-G] 

j=l  i=l 


o 


f 


1 


n . *n  . -1 
1 D 


I [f(i,j)-F]2}2 

j=l  i=l 


o 

g 


n . *n . -1 


n . n . ^ 

I I [g(if j)-G]^}^ 
j=l  i=l 


Speod-of- cloud  motion  and  angle-of -cloud  motion  is  trivially 
calculated  as  shown  in  the  next  section  from  the  displacement 
values  (p ' , q ' ) . 

Computation  of  a cross-correlation  matrix  of  25x25  lag 
values  (p,q)  required  approximately  ten  minutes  on  the  IBM 
System  360  Model  50,  according  to  Leesc  et  al.  [71].  The  com- 
putation time  can  be  greatly  reduced  by  using  fast  Fourier 
transform  methods  to  compute  the  cross-covariance  Cov(p,q). 

Since  the  Fourier  transform  of  the  cross-covariance  of  two 
functions  is  the  product  of  the-  complex  conjugate  of  the  Fourier 
transform  of  one  of  the  functions  with  the  Fourier  transform 
of  tho  other  function,  the  inverse  Fourier  transform  of  this 
product  will  then  yield  the  cross-covariance  function.  Time 


81 


required  to  compute  a 64x64  cross-correlation  matrix  on  the 
IBM  360/50  by  the  fast  Fourier  transform  methods,  search  for 
the  maximum  v^^lue,  and  print  out  the  results  was  only  thirty 
seconds  as  compared  to  ton  minutes  for  the  direct  method. 

Fast  Fourier  transforms  compute  sums  of  products  of  two 
image  matrices  at  different  lag  positions.  These  sums  can  be 
used  to  calculcite  the  dot  product  betv;een  two  vectors  (cross 
correlation)  or  the  distance  botv/ecn  two  vectors  (Euclidean 
norm) . Advantages  of  the  Euclidean  norm  technique  mentioned 
by  Phillips  and  Smith  [72]  are  its  high  degree  of ' flexibility 
and  readily  interprctablc  results. 

The  GOES  simulation  model  for  wind  velocitv  estimation  at 
the  National  Environmental  Satellite  Service  uses  a first- 

guess  displacement  in  combination  with  direct  computation  of 

i 

I 

cross-correlation  coefficients  to  obtain  a cloud-displacement 
vector.  The  computer  program,  which  is  specifically  adapted 
to  operate  as  efficiently  as  fast  Fourier  transform  methods  on 
the  CDC-6600,  calculates  an  initial  set  of  about  100  values  of 
cross-correlation  coefficients  in  the  neighborhood  of  the 
first-guess  location  rather  than  the  entire  set  of  4096  co- 
efficients for  a 64x64  matrix.  If  the  displacement  determined 
from  the  initial  sot  of  100  values  matches  the  first-guess 
field  for  the  cloud-motion  vector  (whicli  is  derived  from  cloud- 
height  information  coupled  with  wind  field  profiles  from  tlie 
National  Meteorological  Center),  the  cloud  displacement  is 
accepted.  The  ini-tial  matching  threshold  was  45  degrees  in 
direction  and  5 meters  per  second  in  speed.  If  the  displace- 
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mont  vector  docs  not  match  the  first-guess  field,  the  entire 
matrix  of  cross-correlation  coefficients  is  computed. 

Although  cross-correlation  could  be  applied  to  the  binary- 
edge  pictures  constructed  by  Bristor's  [72]  binary  matching 
technique  described  in  the  previous  section,  processing  time 
requirements  wore  about  equal  or  slightly  less  for  cloud  dis- 
placement calculation  from  match  scores  of  the  sector  from  the 
first  image  with  every  possible  (or  reasonable)  sector  from 
the  second  imago.  The  origin  of  the  cloud-displacement  vector 
was  defined  as  tlic  center  of  the  first  sector.  The  number  of 
positions  containing  I's  in  both  the  sector  from  the  first 
image  and  a sector  from  the  second  image  formed  the  match 
score.  The  terminus  of  the  cloud-displacement  vector  was  de- 
fined as  the  center  of  the  sector  from  the  second  image  with 
highest  match  score.  The  method  was  critically  dependent 
upon  the  equivalence  of  the  brightness  slice  in  the  tvjo  sec- 
tors. If  it  were  possible  to  calibrate  and  normalize  the  data 
for  solar  illum.ination , methods  such  as  the  binary  matching 
technique  or  ISODATA  technique  which  depend  upon  thresholds  to 
determine  cloud  edges  or  cloud  masses  could  bo  used  operation- 
ally to  compute  cloud-displacement  vectors,  according  to  Leese 
164]. 

A program  called  MOTION,  developed  at  Stanford  Research 

Institute  and  described  by  Endlich  ct  al. [71],  matched  ISODATA 

brightness  centers  in  two  successive  frames  by  iterating  on  a 

fitting  function  F.  F measured  the  distance 

1 

2 2 2 2 

{ (Ax-AXj^)  + (Ay-Ayj^)  + (AB)  } between  the  displacement 
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vector  for  any  two  centers  taken  from  successive  frames  and 


an  input  displacement  vector  (Ax  , Ay,  , 0)  with  x displace- 
mcQt  AXj^ . y displacement  Ayj^  / .and  brightness  displacement  of 
zero.  The  initial  values  AXj^  and  Ay^^  wore  specified  as  the 
median  displacement  in  the  x direction  and  the  median  displace- 
ment in  the  y direction  respectively  taken  over  all  pairs  of 
centers  with  the  first  center  in  picture  1 and  the  second 
center  in  picture  2.  The  pair  of  centers  for  v/hich  F was 
smallest  was  then  considered  to  be  a matched  pair.  The  match- 
ing procedure  continued  for  centers  not  previously  matched  by 
searching  again  for  the  smallest  value  of  F.  The  matching 
procedure  terminated  either  wheii  there  were  no  more  centers  to 
be  matched  or  when  the  value  of  F exceeded  a predetermined 
threshold  value  chosen  in  this  case  to  be  representative  of  a 
departure  from  .the  input  displacement  vector  of  approximately 
lOm-sec  For  the  next  iteration  and  for  every  succeeding 

iteration  of  F values,  the  input  displacement  values  A.Xj.  and 
Ay,  were  chosen  as  average  values  of  x displacements  and  y 
displacemcntr,  for  matched  pairs  only.  In  the  cases  investi- 
gated by  Endlich  et  al.  [71]  , only  three  iterations  were  neces- 
sary for  stable  results  in  pair  matching. 

Approaches  for  selecting  ISODATA  brightness  centers  other 
than  applying  ISODA.TA  techniques  to  the  two  successive  pictures 
independently  were  analyzed.  It  was  found  that  movement  of 
groups  of  clouds  across  boundaries  in  some  cases  presented 
problems  when  using  the  independent  analysis  approach  for 
successive  frames.  It  was  then  decided  to  use  the  centers  on 
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the  I'irr.t  picture  ^is  an  initial  guess  for  centers  on  the  second 
picture;  this  v.'orked  v/ell  if  tlie  general  displacement  was  not 
too  larg''  in  comjiarison  witli  tlie  avcriige  distance  between 
brightness  cent'-rs.  T’ne  approach  finiilly  adopted  used  an 
initial  estimate  of  the  general  displacement  to  reposition  the 
brightness  centers  from  picture  1 as  initial  guess  centers  for 
picture  2.  The  initial  estinic'ito  of  the  general  displcacement 

v.’as  obtained  as  follows.  .Splitting  commands  of  ISODATA  applied 

independently  to  the  two  picture.s  transformed  1 center  into 
2,  then  4,  then  8 centers  for  each  picture.  I’he  MOTION  program 

matched  as  many  as  possible  of  the  8 centers  from  picture  1 

witli  centers  from  picture  2.  The  average  x and  y displacements 
of  these  matched  pairs  determined  the  initial  estimate  of  the 
general  di  s;ilacom.ent . 

Tlio  ISODATA  brightn.ess  centers,  which  represent  cloud 
data  v.’ithin  an  area  having  a certain  radius,  yielded  averaged 
motion  vectors  apprcnria'-.c  for  use  in  synootic-scalo  numer- 
ical f oreCQ.st  ing . Difficulties  v.'ore  encountered  in  aj7plying 
the  ISODATA  technique  to  large,  dense, cloud  bands  such  as  the 
intertroj-ucai  convergence  zone  in  v;hich  a variety  of  cloud 
motions  v.’ith  few  distinct  discernible  features  could  bo  seen. 
The  human  eye  tends  in  those  cases  to  follow  any  individual, 
small,  recognizable  features  which  are  conserved  from  frame 
to  fram.c  rather  tlian  average  the  pattern  into  10-15  brightness 
fields.  To  improve  the  accuracy  of  cloud-displacement  vectors 
for  multi-layered  cloud  regions,  the  addition  of  infrared 
measurements  to  form  a four-dimensional  input  vector  for  tlie 
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ISODATA  progr:iin  was  sucjgcstcd. 

Lo  and  Parikh  [73 J compared  a Fourier  cross-correlation 
method  and  a Fourier  phase  difference  melhod  for  estimation 
of  cloud -d isplaci.'nien t vc'Ctors  on  a variety  of  simuiatod  cloud 
data.  .Simulated  condilionr.  included  clianqe  of  brigliLncss, 
change^  of  si/.e,  cioud  rotation,  c'dgo  effecLs  (cloud  elements 
moved  out  of  tlio  edge  o'"  one  i^icture  but  similar  clouds  did 
not  come  in  at  tfie  ot);.)osito  edge  of  the  second  picture.)  , mi.xcd 
motion,  and  random  'instrumental  noise. 

The  ph,ase  difference  method  is  based  on  the  assumption 
that  the  second  picture  g(x,y)  represents  a linear  shift  of  p 
units  i.n  the  x directio.n  and  q units  in  the  y direction  from 
tlie  first  picture  f(x,y),  i.c., 

y(>:,y)  " f(x-p,y-q)  for  all  x,y 
The  PouriL'i'  transform  H(u,v)  of  the  cross-covariance  function 
h(x,y)  of  tv.’O  pictures  f(x,y)  and  g(x,y)  is  given  by 


li  ( u , V ) = G ( u , V ) * • F ( u , V ) 


where 

G(u,v)*  = com.plc.x  conjugate  of  Fourier  transform  of  the 
sccor.d  picture  g(.x,y) 

F(u,v)  = Fourier  transform  of  the  first  picture  f(x,y) 
From  the  definition  of  the  Fourier  transform  and  the  assump- 
tion that  g(x,y)  = f (x-p,y-q)  , it  can  bo  seen  that 


G(u,v)  = I g(x,y)  e dxdy 


f(x-p,y-q)e  ^ ^ ^ dxdy 


Cr>  CO 


-2ti  i [ (x '+p)  u+ (y  ' +q)  v]  dx'dy' 


= J J f (x ’ ,y ' ) G 

— OO  —00 

If  F{u,v)  is  expressed  in  torms  of  its  amplitude  and  phase 
angle  as 

F(u,v)  = I |F(u,v)  1 le'^'^f 

then 

GMu,v)  = 1 lF(u,v)  1 lc^^2r(pu+qv)-c;,^(u,v)] 

and 

H(u,v)  = G*  (u,v) -F  (u,v)  = 1 ]f  (u,v)  I 1 
The  phase  angle  y^(u,v)  for  H(u,v)  is  given  by 
Vj^(u,v)  = 2:-  (pu+qv) 

The  siiifts  p and  q can  be  determined  from  successive  u and  v 
values  respective] y from  the  relations 
V,^(u  + l,v)  - ,-j^{u,v)  = 2t:P 

Vj^(u,v+1)  - -,q^(u,v)  = 2T!q 

t.he  j'har.e  difference  metl-iod  and  the  cross-correlation 
method  v/erc  applied  to  simulated,  mixed-motion  data  which  vio- 
l^^tcd  file  fiasic  assumption  of  linear  shift 

(i.e.,  g(x,y)  / f(x-p,y-q)  for  any  p,q),  neither  the  phase 
difference  method  nor  the  cross-correlation  method  could  be 
analysed  to  obtain  true  estimates  for  both  types  of  motion. 
Filtering  in  the  frequency  domain  or  thrcsliolding  in  the  spatial 
domain  were  suqgc.stcd  as  possible  approaches  for  future  research. 
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3.4.3.  Automatic  Computation  of  Wind  Velocity 


Tlie  output  from  automatic  techniques  for  tracking  cloud 

target.n  is  a displcicement  vector  (p',q').  Speed  |c]  cind 

direction  0 of  cloud  motion  are  calculatc.'d  as  follows: 

1 

|c|  = [ ([^'  -Ax)^  -t-  (.j'  -Ay)^]^ 

At 


0 


arctan 


p ' -Ax 
q ' • Ay 


where 


Ax,  Ay  = .sa'iipl  i ncj  intervals  of  the  input  picture  matrix 
in  tli(.'  ea.s  (. -v.-est  and  nor  tli-.sou  t h directions, 
rLuspecti  vc'J  y 

At  = t i ::ie  interval  between  th.c  two  jjictures 

The.'  automatr'd  mofk.ds  of  tlic'  At'itional  bnvi.ronmental 
Sc'itellit(^  Service  (NE.SS)  for  estimation  of  low-J.cvel  cloud 
motion  ve'ctoiT.  h:ive  incortx.n-atcd  an  automatic  ’p,ost-cditor  as 
a jnajor  part  of  llie  quality  control  unit..  T)ie  post-editor  of 
the  ACI.OWi  r'.odei  of  h'lASS  rejects  a low-level  cloud  motion 

\ectv'r  .ir.  inva.l  id  if  il  sa.t  isflc's  either  of  tlie  following 
climatological  cri  tori. a : 

(1)  it  5'.  r.K'uin.i  Lude  is  greater  t.han  40  kts, 

(2)  it  has  a westerly  componcnit  and  is  located  between 
20"N  and  2 0'’S. 

The  post-editor  cif  the  GOES  simulation  model  of  NESS  considers 
criteria  such  as  displacement  from  tlie  first-guess  field,  hori 
zontal  consistency  between  neighboring  cloud-motion  vectors, 
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and  gradient  around  the  primary  maxima  of  the  cross-correlation 
field  for  determination  of  the  validity  of  lov;-levol  cloud- 
motion  vectors. 
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. Conclusions  and  RccontTiendations 


The  purpose  of  this  study  was  to  assess  the  current  status  of  pattern  rec- 
ognition techniques  for  severe  storm  features  and  to  consider  the  potential 
applicability  of  pattern  recognition  techniques  to  the  automatic  identifica- 
tion and  prediction  of  severe  weather  phenomena  from  satellite  data.  In  the 
previous  sections  of  this  report,  wo  have  reported  in  some  detail  on  the  na- 
ture of  severe  storm  features  obtainable  from  visible  and  infrared  satellite 
data.  We  have  also  described  the  past  develonments  and  current  status  of 
automatic  techniques  related  to  cloud  classification  and  weather  phenomena. 

In  this  section,  we  summarize  our  conclusions  and  recop'.mcnd  directions  which 
appear  -fruitful  for  further  investigation. 

4.1.  Conclusions 

In  the  last  decade  automatic  pattern  recognition  and  image  processing  tech- 
niques for  analysis  of  meteorological  satellite  data  have  been  directed  pri- 
marily toward  solution  of  cloud  classification  and  wind-velocity  estimation 
problems.  The  resolution  and  spectral  characteri sti cs  of  the  meteorological 
satellite  data  controlled  the  type  of  information  which  could  be  automatically 
extracted  from  the  data  and  the  nature  of  apnli cable  pattern-recognition  tech- 
niques. For  coarse  resolution  data,  point-by-point  classification  schemes  were 
developed.  With  the  availability  of  high-resolution  scanning  radiometer  data, 
decision-tree  models  based  on  extraction  of  cloud-type  features  wore  constructed. 
Critical  thresholds  for  cither  poi nt-by-poi nt  classification  schemes  or 
decision-tree  models  were  determined  either  from  design  samples,  from  theoreti- 
cal considerations,  or  from  conventional  meteorological  sources.  Cloud-shape 
and  texture  features  could  be  determined  from  visible  data  only.  Successful 
cloud-type  classification  schemes  required  multispoctral  satellite  data. 
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Geosynchronous  satellite  data  formed  the  orimary  input  to  automatic  wind- 
velocity  estimation  programs.  Cross-correlation  techniques  were  applied  to 
two  successive  windov/s  of  either  infrared  or  visible  data  to  determine  cloud 


displacement  vectors.  Techniques  to  extract  wind  estimates  from  displace- 
ments confined  to  individual  thermal  slices  are  currently  under  investiga- 
tion at  NES';. 

The  automatic  cloud-classification  models  and  automatic  wind-velocity  esti- 
mation models  have  been  extensively  tested  only  over  oceanic  regions.  The  non- 
homogeneity of  the  land  surface  as  a background  coupled  with  an  increase  in  the 
influence  of  local  meteorological  phenomenon  which  result  in  non-advective  cloud 
motions,  complicates  the  problem  of  automatic-cloud  type  and  wind-velocity  analysis 
over  land  areas.  The  success  of  these  systems  on  a global  scale  has  yet  to  be  tested. 

Automatic  severe  weather  identification  and  forecasting  models  place  even 
more  stringent  requirements  than  automatic-cloud  type  or  wind-velocity  models 
on  the  amount  and  nature  of  meteorological  satellite  data  which  must  be  con- 
sidered by  the  system.  In  assessing  the  potential  of  pattern-recognition 
techniques  to  identify  severe  storm  features  from  satellite  data,  we  observed 
that  severe-storm  features  could  generally  be  grouped  into  three  major  cate- 
gories : 

1)  wi nd-depondont  features 

2)  enhanced-! nf rared  features 

3)  visual  features. 

Specific  severe  weather  features  in  each  of  the  above  categories  are  listed  in 
Table  Wind-dependent  features  and  enhanced-infrared  features  (with  the 
possible  exception  of  steepest  temperature  gradient  and  top  heights)  require 
time  sequences  of  data  to  track  the  evolution  of  the  feature. 
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Table  2.  Categories  of  Severe-Storm  Features. 


Relationships  of  specific  severe-storm  features  to  either  wind  data,  en- 
hanced infrared  data,  or  visible  data  suggested  by  Table  1 represent  one  and 
only  one  approach  to  recognition  and  quantitative  measurement  of  severe-storm 
features.  Alternative  relationships  exist.  Automatic  techniques  to  identify 
pendant-sfiaped  colls,  for  example,  could  be  formulated  based  on  cell-shape 
measurements  rather  than  on  veer  of  wind  with  height.  Although  overshooting 
tops  were  categorized  under  enhanced  infrared  features,  the  resolution  of  in- 
frared SMS  data  is  too  coarse  to  delineate  individual  overshooting  tops.  Out- 
lines of  individual  tops  are  best  observed  in  higher  resolution  visible  imagery 
for  a low-sun  angle.  The  major  disadvantage  of  using  visible  data  to  track 
overshooting  tops  is  that  a relatively  expensive  normalization  procedure  for 
solar  zenith  angle  would  be  necessary  before  statistical  relationships  between 
top  brightness  and  severe-weather  phenomena  cculd  be  develoned. 

Our  general  conclusion  in  examining  the  current  status  of  severe-storm  iden- 
tification is  ttiat  it  is  an  area  of  individual,  subjective  experience  which 
exists  in  the  form  of  verbal  descriptions  of  various  features.  As  a result, 
thene  arc  differences  of  opinion  among  meteorologists  concerning  the  nature 
of  tlie  data  necessary  for  identification  of  particular  severe  weather  fea- 
tures and  the  role  of  automatic  processing  in  forecasting  severe  weather  phe- 
nomena. Although  most  of  the  responses  to  our  severe  weather  questionnaire 
indicated  a need  for  automatic  processing  of  high  resolution,  time  sequences 
of  visible  and  infrared  meteorological  satellite  data  in  an  interactive  en- 
vironment, the  blend  of  man-machine  interaction  envisioned  within  these  sys- 
tems varied. 

On  the  one  hand,  it  is  said  that  automatic  procedures  cannot  hope  to  com- 
pote with  computerized  methods  whO'"  dealing  with  local  phenomena  over  small 
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regions.  A meteorologist  can  cuickly  assess  ttie  relevance  of  ground-based 
reports,  teletype  reports,  radar,  and  satellite  features  and  bring  to  bear 
his  expertise  on  the  limited  area  of  concern.  In  this  view,  the  role  for  au- 
tomated techniques  is  in  dealing  wifh  more  global  situations  where  vast  amount 
of  data  are  to  be  handled  in  order  to  define  the  small  percentage  of  areas  on 
which  a human  should  concentrate. 

On  the  other  hand,  image  processing  of  large  amounts  of  satellite  data  on  a 
sequential  computer  is  not  economically  feasible  at  the  present  time.  The 
general  procedure  to  handle  large  amounts  of  data  is  to  partition  the  data 
into  smaller  regions  (windows)  and  base  decisions  on  local  operations  per- 
formed on  each  window.  This  procedure  fails  to  give  satisfactory  results 
v/hen  contextual  information  from  nearby  windows  is  necessary  for  a decision 
or  when  given  features  lie  in  two  different  windows.  As  a consequence,  many 
meteorologists  prefer  to  use  automatic  techniques  to  exhaustively  examine 
small  regions  which  they  have  selected  in  order  to  detect  and  measure  subtle 
differences  and  changes  which  are  not  easily  measured  by  the  human  visual 
system. 

It  would  seem  that  a more  careful  study  of  the  interaction  between  human 
and  machine  processing  of  weather  data  will  be  facilitated  by  whatever  attempts 
are  made  in  automating  some  of  the  severe  storm  features  in  Table  As  a 
first  step,  descriptions  of  severe  weather  features  should  be  standardized 
so  as  to  allow  automatic  extraction  whenever  feasible.  For  wind-dependent 
features  in  Table  algorithms  to  extract  low-level  and  upper-level  wind 
vectors  over  oceanic  areas  and  to  determine  divergence  and  vorticity  (Thomasell 
[1977])  arc  being  investigated  at  KESS.  Our  extensive  search  of  the  litera- 
ture and  our  questioning  of  individuals  working  on  different  aspects  of  me- 
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teorological  data  processing  revealed  no  existing  algorithms  for  automatic 
extraction  of  enhanced  infrared  and  visible  imagery  of  severe  storm  features. 
While  many  pattern  recognition  algorithms  and  image  processing  algorithms 
exist  for  the  various  lov;-level  operations  that  would  be  involved  in  extrac- 
tion of  these  features  (see  Section  3),  a comparative  study  of  the  applica- 
bility of  relevant  algorithms  does  not  exist.  To  determine  the  best  algo- 
rithms relative  to  these  features  would  require  experimentation  with  real  data. 

As  mentioned  earlier,  in  the  pattern  recoonition  and  image  processing  litera- 
ture almost  no  attention  has  been  paid  to  the  extraction  of  features  which 
evolve  with  time.  This  time  variation  of  enhanced  infrared  severe  storm  fea- 
tures has  not  been  attacked  previously  using  pattei  r.  recognition  techniques, 
and  any  demonstration  of  feasibility  on  these  features  would  lend  credence 
to  pattern  recognition  approaches  for  severe-storm  pattern  recognition.  If 
these  features  could  be  automatically  observed  and  extracted,  this  would  sub- 
stantially aid  in  rapid  decision  making  concerning  severe  storm  activity  and 
testing  correlations  between  measurements  from  satellite-derived  severe-storm 
features  and  from  ground-based  reports  on  severe  weather  phenomena.  For 
this  reason,  it  is  worthwhile  to  consider  the  development  of  automatic  tech- 
niques . 
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4.2.  Reconifiicnda  ti  on  s 

It  is  recommended  that  an  attempt  be  made  to  quantify  current  verbal  de- 
scriptions of  severe-storm  features  in  order  to  minimize  ambiguity  of  defini- 
tion atid  to  facilitate  (a)  comoarativc  studies  of  the  correlation  of  the 
various  features  to  severe-weather  phenomena  and  (b)  the  automatic  extraction 
of  those  features  whicli  correlate  well  with  specific  severe-weather  phenomena. 
Implementation  of  this  recommendation  will  require  close  and  extended  inter- 
action with  meteorologists  active  in  severe-storm  identification  and  predic- 
tion from  satellite  data. 

The  key  to  the  automation  of  many  of  the  potentially  useful  severe-storm 
features  lies  in  being  able  to  detect  and  recognize  patterns  which  evolve  in 
time.  It  is  imoortant  to  establish  os  early  as  possible  the  feasibility 
of  automatically  recognizing  such  patterns.  It  appears  that  it  may  be  pos- 
sible to  recognize  some  of  these  patterns  using  limited  scan  infrared  data. 

In  the  event  that  it  should  Drove  feasible  to  automatically  extract  time- 
varying  patterns  as  recommended  above,  it  would  be  desirable  to  conduct  a 
comparative  study  of  curve-detection  operators  and  line-detection  operators 
for  identification  of  mesoliigh-boundary  interactions  and  low-level  cumulus 
cloud-line  feeder  bands.  In  lino  with  the  above,  it  is  recommended  that  ini- 
tially work  bo  done  on  merging  or  splitting  cells  in  infrared  data  enhanced 
in  a manner  similar  lo  that  currently  being  used  for  human  recognition.  For 
this  investigation,  data  which  contains  well-defined  events  of  the  type  to  be 
recognized  will  be  needed. 

Algorithms  for  wind-velocity  estimation  and  calculation  of  divergence  and 
vorticity  fields  over  ocean  areas  are  being  develoned  at  NESS.  There  are  cur- 
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rently  differences  of  opinion  among  meteorologists  on  the  feasibility  of  using 
automatic  techniques  to  calculate  cloud  displacements  which  represent  the  am- 
bient wind  flow.  The  applicability  of  automatic  algorithms  to  determine  wind- 
veclocity  vectors  from  cloud  displacements  over  land  has  not  been  determined. 

In  the  event  that  they  are  applicable,  these  programs  could  be  incorporated 
into  automatic  metliods  to  extract  wind-dependent  severe-storm  features. 

Decision  trees  (similar  to  that  of  Scofield's)  which  make  use  of  a number  of 
features  to  determine  specific  severe-storm  activity  should  be  developed.  In 
developing  these'  trees,  particular  attention  should  be  given  to  features  whose 
extraction  is  likely  to  be  successfully  automated.  Interactive  programs  for 
the  design  of  such  decision  trees  should  also  be  considered. 
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APPENDIX 


Severe  Storm  Questionnaire 


OUESTIONNAIRi: 


Please  send  responses  to: 

Professor  Laveen  Kanal 
Laboratory  for  Pattern  Analysis 
Departnent  of  Computer  Science 
University  of  Maryland 
College  P;irk,  Maryland  20742. 


YOUR  NA-'lE:  

ADDRESS  & TEL.  NO. 


1.  What  are  the  cc,*'C^oj. ies  of  severe  weather?  (Some  examples 
are  listed  below.  What  additional  categories  should  be 
listed? ) 

1.  Severe  windstorm 

2.  Heavy  r?;.instorm 

3.  Hailstorm 

4 . Tornado 

5.  

6.  

7.  

S. 

9. 


2.  Correspondina  to  the  above  categories,  please  list  the  im- 
portant characteristics  (features)  of  each  of  those  categories 
which  you  think  may  be  obtainable  from  SMS  satellite  data;  e.g., 
splitting  colls,  pendant  shaped  cells,  overshooting  thunder- 
storm tops,  right  moving  storms,  low-level  cumulus  feeder 
bands,  mesohigh  boundary  interactions,  steepest  temperature 
gradient,  speed  of  coll  motion,  vertical  rate  of  cell  growth. 
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QUEST  TOtiU  AIRE 
Payo  2 


2.  (Continued) 

pattern  of  upper  level  and  lower  level  convergence  and  di- 
vergence. If  possible,  please  list  the  features  for  each 
category  in  the  order  of  their  importance  for  characteriz- 
ing each  category. 


3.  Which  of  the  above  characteristics  do  you  feel  a’"c  amenable 
to  computer  e.':traction  from  visible  and  IR  satel.ite  data? 


QUESTIONNAIRE 
1?  ag  a J 


4.  For  each  of  the  charac'.  or istics  listed  under  question  3,  wo 
would  like  to  l;nov.’  your  opinion  as  to  whether  or  not  that 
characteristic  can  bo  extracted  from  a single  IR  image.  If 
not,  arc  visible  imagery  and/or  a sequence  of  visible  and/ 
or  IR  images  required?  If  so,  how  many  images  and  at  what 
intervals  do  you  suggest? 


5.  Do  you  know  of  algorithms  which  exist  or  are  under  develop- 
ment for  computer  extraction  of  such  features?  If  so,  please 
provide  references  to  the  algorithms  and  to  the  laboratories 
or  persons  who  can  be  contacted  concerning  them. 


QUESTIONNAIRE 
Page  4 


6.  Arc  there  algorithms  which  you  feel  should  bo  implemented 
either  for  preprocessing  meteorological  satellite  data  or 
for  analyzing  specific  severe  storm  areas?  Do  you  feel 
that  such  areas  will  have  to  be  pin  pointed  by  a m.eteorolo- 
gist  using  a man-machine  interactive  system? 


7.  Do  you  have  specific  suggestions  on  how  computers  can  best 
help  in  severe  storm  weather  pattern  recognition  using  sa- 
tellite data? 
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